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Abstract 

We study the chromatic polynomial Pg{q) foi' rn x n square- and triangular-lattice 
strips of widths 2 < m < 8 with cyclic boundary conditions. This polynomial gives the 
zero-temperature limit of the partition function for the antiferromagnetic g-state Potts 
model defined on the lattice G. We show how to construct the transfer matrix in the 
Fortuin-Kasteleyn representation for such lattices and obtain the accumulation sets of 
chromatic zeros in the complex g-plane in the limit n — > 00. We find that the different 
phases that appear in this model can be characterized by a topological parameter. We 
also compute the bulk and surface free energies and the central charge. 



Key Words: Chromatic polynomial; antiferromagnetic Potts model; triangular lattice 
square lattice; transfer matrix; Fortuin-Kasteleyn representation; Beraha numbers; confor 
mal field theory. 



1 Introduction 



The two-dimensional Potts model conceals a very rich physics, in particular in the anti- 
ferromagnetic regime. Let G = {V,E) be a finite undirected graph with vertex set V and 
edge set E. The gr-state Potts model is initially defined for q a positive integer in terms of 
spins ax — 1,2, ... ,q living on the vertices x &V. Its Hamiltonian reads 

^(M) = -J Yl ^K'^.)' (1-1) 

(xy)£E 

where 6 is the Kroncckcr delta. Fortuin and Kasteleyn [1] have shown that the partition 
function at inverse temperature /? can be rewritten as 

ZgM = J^g'^^^VI, (1.2) 

E'CE 

with V — e^"^ — 1. The sum runs over the 2l^l subsets E' C E, with k{E') being the number of 
connected components (including isolated vertices) in the subgraph {V, E'). We now promote 
(1.2) to the definition of the model, which permits us to consider q as an arbitrary complex 
number. 

We shall here be interested in the zero-temperature antiferromagnetic case J = — oo (i.e., 
V — —1), whose restriction to g a positive integer can be interpreted as a coloring problem. 
In this case, Pg{q) = Za{q,—1) is known as the chromatic polynomial. Its evaluation for 
a general graph is a hard problem in the sense that the best algorithms currently available 
have a worst-case execution time which grows exponentially with the number of sites of 
the graph. For strip graphs (or, more generally, recursive families of graphs) the chromatic 
polynomial Pg„ (?) of any member of the family G„ can be computed from a transfer matrix 
T and certain boundary condition vectors u and v: 

Pg„ = u^-T'^-v. (1.3) 

Thus, the computation time grows as a polynomial in n; it does however still grow expo- 
nentially in the strip width (due to the L dependence of dimT). Thus, it is necessary to 
introduce new algorithms to be able to handle as large L as possible. For example, identify- 
ing symmetries may give an equivalent transfer matrix of lower dimension (in this respect, 
the Fortuin-Kasteleyn representation makes manifest the Sq symmetry of the original spin 
Hamiltonian) . 

It is important to remark that the Potts spin model has a probabilistic interpretation 
(i.e., has non- negative Boltzmann weights) only when g is a positive integer and f > — 1. The 
Fortuin-Kasteleyn random-cluster model (1.2), which extends the Potts model to non-integer 
q, has non-negative weights only when q >0 and v >0. In all other cases, the model belongs 
to the "unphysical" regime (i.e., some weights are negative or complex), and some familiar 
properties of statistical mechanics need not hold. For example, even for integer q > 2 and 
real v in the range —l<v<0, where the spin representation exists and has non- negative 
weights, the dominant eigenvalue Ai of the transfer matrix in the cluster representation need 
not be simple (i.e., the eigenvector may not be unique); and even if simple, it may not 
play any role in determining Za^{q,v) because the corresponding amplitude may vanish. 
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Both these behaviors are of course impossible for any transfer matrix with positive weights, 
by virtue of the Perron-Frobenius theorem. It is important to note that the dominant 
eigenvalues in the cluster and spin representations need not to be equal, because one of them 
may have a vanishing amplitude (see below). In any case, the Potts model can thus make 
probabilistic sense in the cluster representation at parameter values where it fails to make 
probabilistic sense in the spin representation, and vice versa. It is worth mentioning that for 
q = 4cos^(7r/n) with integer n = 3,4, . . . and planar graphs G, the Potts-model partition 
function Zdq, v) admits a third representation, in terms of a restricted solid-on-solid (RSOS) 
model [2, 3]. 

But even in the most general case (i.e., with arbitrary real, or even complex q and v) we 
consider the study of the random-cluster model (1.2) fully legitimate. Note in this context 
that there is a long story of studying statistical-mechanical models in "unphysical" regimes 
(e.g., the hard-core lattice gas at its negative-fugacity critical point [4, 5, 6, 7, 8, 9, 10, 11, 
12, 13] or the Yang-Lee edge singularity [14, 15, 16, 17, 18, 19, 20]). In particular, Baxter 
[21] studied the zero-temperature limit of the triangular-lattice Potts antiferromagnet, which 
belongs to this "unphysical" regime. Generally speaking, the "unphysical" regime can be 
understood using the standard tools of statistical mechanics with appropriate modifications. 
In particular, although a rigorous proof is still lacking, conformal field theory (CFT) seems 
to apply also in the "unphysical" regions (sec e.g., [19, 20] for the CFT description of the 
Lee-Yand edge singularity). In general, "unphysical" regions are described by non-unitary 
field theories. 

It is well established that antiferromagnetic models in general, and the chromatic polyno- 
mial in particular, are very sensitive to the choice of boundary conditions. Indeed, different 
choices may lead to quite different thermodynamic limits. Following our earlier publications 
[22, 23, 24], where we studied free or periodic boundary conditions in the transverse (space- 
like) direction and free boundary conditions in the longitudinal (time-like) direction, we here 
subject the chromatic polynomial to cyclic boundary conditions, i.e., free in the transverse 
direction and periodic in the longitudinal direction. For a transfer matrix written in the Potts 
spin basis, cyclic boundary conditions can be easily implemented by setting Pq^ = TrT". 
The situation is, however, less straightforward in the Fortuin-Kasteleyn representation, and 
we are obliged to employ a transfer matrix that is different from — and more complicated 
than — the ones used in [22, 23, 24]. 

Our study is of interest for understanding the general properties of the so-called Berker- 
Kadanoff (BK) phase in the antiferromagnetic Potts model [26]. Generically, the two- 
dimensional Potts model on a given regular lattice possesses a curve vfm{q) > of fer- 
romagnetic phase transitions which are second-order in the range < g < 4 [27]. Along 
this curve the thermal operator (which corresponds to taking v away from its critical value 
vpuiq)) is relevant and it becomes marginal in the limit g ^ of spanning trees. The 
analytic continuation of the curve VpuiQ) into the antiferromagnetic regime yields a second 
critical curve vbk{q) < with < g < 4 along which the thermal operator is irrelevant. 
Therefore, for each q e (0,4), the critical point i'bk(?) acts as the renormalization-group 
attractor of a finite range of v values: this is the BK phase. 

The above scenario has been verified in detail for several regular lattices. In particular, 
for the square and triangular lattices the corresponding models are Bethe Ansatz soluble 
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[28]^ and the analytic forms of i'fm(?) and i'bk(?) are exactly known. The scenario is also 
inherent in the Coulomb gas approach to critical phenomena [29] , which, albeit less rigorous, 
confirms that the general properties of the BK phase are lattice independent. 

The best understood case is that of the square lattice. Here, vpm{q) = and 
VBK{q) = [28, 26]. There are also two mutually dual (and hence equivalent) curves 

of antiferromagnetic transition points [30], v±{q) — —2 ± •^^4 — q for q e [0,4], which are 
believed to form the boundaries of the BK phase [26]. Exactly on v±{q), the model exhibits 
a different type of critical behavior [26, 31]. 

More precisely, for a given width L and q G I{L) belonging to a suitable interval I{L), 
there exists mutually dual couplings v±{q,L) so that for V-{q,L) < v < L), the 

dominant eigenvalue XBKiv) in the cluster representation is simple. As L — > oo, one has 
v±{q,L) — ^ v±{q) and limi^oo ^iL) contains [0,4]. The finite-size scaling of \bk{v), and of 
its adjacent scaling levels, determines the critical exponents of the BK phase. However, for 
certain special values of q (the generalized Beraha numbers, see below) the amplitudes of 
Xbk{v), and of many subdominant eigenvalues, vanish; yet other eigenvalues are degenerate 
in pairs but have amplitudes of opposite sign and thus, cancel out from the partition function. 
The remaining eigenvalues correspond to those of the relevant RSOS representation (or spin 
representation, when g is a positive integer). On the other hand, for v ^ [v^{q, L),v^{q, L)] 
the dominant eigenvalue in the cluster representation is again simple, but is a different 
analytical expression of v, and its amplitude does not vanish at the special values of q. 

On the triangular lattice, i'fm(?) is the uppermost branch {v > 0) of the curve + Sv'^ — 
q — [32] . We know of numerical evidence (to be published elsewhere) that the middle and 
lower branches of this same curve may be respectively VBK{q) and the lower boundary f _(g) of 
the BK phase. Less is known about the upper boundary v+{q) of the BK phase. If we assume, 
in analogy with the square-lattice case, that it is the first locus of partition function zeros 
encountered in the antiferromagnetic regime upon decreasing the temperature from infinity, 
one has ^+(0) = 0, and the slope ■^v+{q) ~ —0.1753 ± 0.0002 is known numerically 

[33]. Moreover, another Bethe Ansatz soluble case [21], v — —1, provides evidence that 
'^+(Q'o) = ~1) where go will be defined below. 

In the present publication we study the chromatic polynomial with cyclic boundary con- 
ditions for the square and the triangular lattices. In both cases the chromatic line v = —1 
passes through a finite part of the BK phase, whence our results will enable us to examine the 
role of cyclic boundary conditions for this phase. We shall obtain strong evidence that with 
cyclic boundary conditions the BK phase splits into an infinite number of distinct phases, 
each of which has a precise topological characterization in terms of the Fortuin-Kasteleyn 
clusters. However, the chromatic line only intersects a finite number of these phases. For 
the square (resp. triangular) lattice chromatic polynomial we thus find a total of three (resp. 
seven) phases. 

In particular, we examine the loci of zeros of PoLxni^i strips of width L and length n, 

^Baxter [28, 30] has shown that the two-dimensional Potts model is equivalent to a homogeneous six- 
vertex model on certain curves of the {q,v) plane, including vpM{q) and wbk(9) for both the square and 
triangular lattices. This homogeneous six- vertex model is soluble using the Bethe Ansatz [28]. In addition, 
Baxter also showed that the square-lattice model on the mutually dual curves v± (q) is equivalent to a certain 
class of inhomogeneous six- vertex models that can nevertheless be solved by an extension of the Bethe Ansatz 
[30]. 
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in the plane of complex q. When n — > oo, the accumulation points of these zeros form either 
isolated limiting points (when the amplitude of the dominant eigenvalue vanishes) or con- 
tinuous limiting curves Bl (when two or more dominant eigenvalues become equimodular). 
These latter curves constitute the boundaries between the different phases just described. 
We refer to [22, 23, 24] for details, and for the relation to the Beraha-Kahane-Weiss theo- 
rem [34] . The rather unusual phase diagram referred to above is hnked to a difference of the 
L ^ oo limits of Bl with free/cyhndrical [22, 23, 24] and cyclic boundary conditions (this 
work) . 

The curves Bl will in general intersect the real g-axis in a number of values, and we define 
qo{L) to be the largest of these values. Also, let qq — Huil^oo Qo{L) . Based on our earher 
studies [22, 23, 24] we expect for a given lattice to be independent of (reasonable) boundary 
conditions. The physical scenario outlined above would then imply that wc should observe 
the BK phase for < q < qo{L). We expect go(sq) = 3 for the square lattice (corresponding 
to ^+(3) = —1). For the triangular lattice, Baxter [21] originally predicted qo{tn) pa 3.81967. 
We have recently revisited the analysis of Baxter's equations [24], arguing in particular that 
his prediction for go(tri) is erroneous, the correct value being go(tri) —2 + ^/3 3.73205. 

The results available in the literature concerning square- and triangular-lattice strips with 
cyclic boundary conditions are limited to widths L < 5. The results for the square lattice 
were obtained for L = 2 by Biggs, Demerell and Sands [35], for L = 3 by Shrock and Tsai 
[36], and for L = 4, 5 by Chang and Shrock [37, 38]. The results for the triangular lattice 
were obtained for L = 2 by Shrock and Tsai [36], and for L = 3, 4 by Chang and Shrock 
[39]. When writing up this paper we learned that Chang and Shrock have extended their 
results to the case L = 5 [41]. 

In this paper we present new exact results for the chromatic polynomial of strip graphs 
of widths L = 6. We have also computed independently the exact formula for the triangular- 
lattice strip graphs of width L — 5; furthermore, we present the computation of the limiting 
curve for these strips.^ For L = 7, 8 we have exactly computed some blocks of the transfer 
matrix from which all the eigenvalues can be inferred (modulo some conjectures that are 
very plausible though unproven, see Properties 1 and 3 in Section 2.3). Unfortunately some 
of these blocks are so large than we have been unable to get the analytical expressions for the 
eigenvalues, as we have done for L < 6. From the exact symbolic expression for these blocks, 
we can nevertheless obtain essentially-exact relevant physical information, such as numerical 
values for the eigenvalues and the limiting curve B. The only source of error in these results 
is the limited (but very high) machine precision of the numerical computations performed 
with MATHEMATICA. We have taken the symbolic computations reported in this paper as 
far as possible, using several months of CPU time. However, because of the exponential 
character (in the width L) of this problem, going beyond L — 8 will require an extremely 
large amount of CPU time and RAM memory. 

Our approach is based on a transfer-matrix formulation of the Fortuin-Kasteleyn repre- 
sentation, which is equivalent to the methods of Shrock and collaborators (See e.g. [40] and 
references therein).^ Large parts of our analysis has been automatized and carried out by 

^The interested reader can find th.e exact expression of the eigenvalues for L < 6 in the MATHEMATICA 
files available as part of the electronic version of this paper in the cond-mat archive. 

^The two approaches are however formulated in somewhat different terms, so the equivalence is only made 
clear by a detailed comparison. 
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a computer, thus allowing to access larger widths L. On the other hand, in our approach 
wc find a natural candidate for an essentially "conserved quantum number" : the number of 
bridges^ of a given connectivity state (see Section 2.2 for details). This number is the link 
to the quantum-group theory worked out by Pasquier and Saleur [25, 26]. Finally, let us 
mention a different approach developed by Noy and Ribo [42] and also based on a transfer- 
matrix formalism. They handle cyclic boundary conditions by introducing a new type of 
operation: edge deletion. 

The outline of the paper is as follows. In Section 2 we present our method for obtaining 
the transfer matrix for a square-lattice strip with cyclic boundary conditions. We describe 
the structural properties of the transfer matrix and its block structure. Section 2.3 also 
provides an alternative (and simplified) method for constructing the transfer matrix which 
has proved very useful for large-width computations. In Section 3 we present the results 
for square-lattice strips up to widths L = 8. In the following Section we adopt our method 
to triangular-lattice strips with cyclic boundary conditions, and the results are given in 
Section 5 for widths up to L = 8. Finally, in Section 6 we discuss the phase diagram of the 
square- and triangular-lattice Potts antiferromagnet at zero temperature. We compute the 
free energy and the central charge as a function of q. Finally, in Appendix A we show an 
alternative way of computing the transfer matrix for a cyclic square-lattice strip. 

2 Transfer matrix for square- lattice strips with cyclic 
boundary conditions 

We want to build the transfer matrix for a square-lattice strip of width m and arbitrary 
length n with cyclic boundary conditions, i.e., free along the transverse direction and periodic 
along the longitudinal direction (Figure 1 shows an example of size 4x2). As in previous 
work [22, 23, 24], we would like to express the transfer matrix T as a product of two matrices 
H and V, representing respectively the horizontal and vertical bonds of the lattice 



For each lattice strip width m we have a different transfer matrix. When free boundary 
conditions are used in the longitudinal direction (i.e., for a strip with free or cylindrical 
b.c), the transfer matrix in the Fortuin-Kasteleyn representation has been explained in 
detail in [22, Section 3]. The basic idea is that the state of the system is characterized by 
a connectivity pattern {vp} of the sites {!,..., m} of the top row, and the transfer matrix 
T(m) acts by attaching one new row above the current top row. 

By contrast, when periodic boundary conditions are used in the longitudinal direction, 
it is insufficient to keep track of the connectivity of the top row alone, as we will eventually 
need to identify the top and bottom rows. Rather, as first sketched in [22, Section 7.4], the 
idea is to keep track of the connectivity pattern of the 2m sites constituted by the current 
top row and the bottom row. Let us call these sites 1, 2, . . . , m and 1', 2', . . . , m', respectively 
(see Figure 1 for an example with m = 4). Initially the top and bottom rows are identical. 
We then enlarge the lattice by adding successive new rows to the top of the lattice, exactly as 

^The number of bridges is equivalent to what Refs. [41, 40] call the level. 
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in [22, Section 3.2]; the join and detach operations act on the sites of the top row, with those 
of the bottom row simply "going along for the ride" . At the end, when we have obtained a 
lattice with n + 1 rows, we identify the top and bottom rows. 

Let us now explain this formalism in greater detail. We characterize the state of the 
top and bottom rows by a connectivity pattern {vp}, which is associated to a partition 
V — {Pi, . . . , Pk} of the set {1, . . . , m, 1', . ■ ■ , "t.'}. It is convenient to introduce the operators 

= vl + D, (2.2) 

Qx,y = I + vS^y (2.3) 

acting on {v-p} (see [22] for more details). Here, / is the identity; D^, is a detach operator 
that detaches site x from the block it currently belongs to, multiplying by a factor g if x is 
currently a singleton; and S^y is a join operator that amalgamates the blocks containing x 
and y, if they were not already in the same block. 

The matrices H(m) and V(m) then act on the top row as follows: 

m— 1 

H(m) = n^M+i (2-4a) 

1=1 

m 

V(m) = (2.4b) 

1=1 

We can similarly define matrices H'(to) and V'(m) that act on the bottom row. 

The chromatic polynomial of a square-lattice strip of width m and length n with cyclic 
boundary conditions is given by 

Pmxn(g) = ■ H(m) • l{mY ■ Via (2.5) 

The left u and right Vid vectors differ from those found for free and cylindrical boundary 
conditions. The right vector Vid now denotes the partition 

= {{l,l'},{2,2'},...,{m,m'}} (2.6) 

(i.e., we start with the top and bottom rims identified). The left vector in (2.5) acts on a 
connectivity state v-p as follows 

u^.vr ^ gl^'l (2.7) 

where \V'\ denotes the number of blocks in the connectivity v-p' obtained from vp by the 
action of several join operators J-^^j,: 

Vp' = (^fj-'^'^'j '^^ (2.8) 

In other words, acts on Vp by identifying the top and bottom rows, assigning a factor of 
q to each block in the resulting partition. 

In this paper we are mainly concerned by the chromatic-polynomial case {v — —1). In 
this case, the matrix H is a projector (i.e., = H). Thus, we can use instead of the transfer 
matrix T [c.f., (2.1)] and the usual connectivity states Vp, the modified transfer matrix 

T(m) = H(m) ■ V(m) • H(m) (2.9) 
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and the basis vectors 

w-p = H(m) • v-p (2.10) 

Note that w-p — iiV has any pair of nearest neighbors of the top row in the same block. 
As we start with the state Wid, and at the end we are identifying the top and bottom rows, 
it is useful to work directly on the modified connectivity basis 

w-p = H'(m) • H(m) • v-p (2.11) 

This choice implies that wp = also if V contains any pair of nearest neighbors of the 
bottom row in the same block. For simplicity, we hereafter drop the tilde to denote the 
modified transfer matrix (2.9), and the hat to denote the modified connectivity basis (2.11). 
In this case, the chromatic polynomial can be written as 

Pm.nil) = ■ T(m)" • W,d (2.12) 



2.1 Structural properties of the chromatic polynomial 

The first observation is that the square lattice with cyclic boundary conditions is a 
planar graph: thus, only non-crossing partitions P of {1, . . . , m, 1', . • • , can occur. As 
we have transverse free boundary conditions, the number of such connectivities is given by 
the Catalan number 

- ^iM) 

Their asymptotic behavior is given by 

= 16"^m-=^/^(87r)-V2[i + 0(i/^)] (2.14) 

These numbers are listed in the second column of Table 1 up to m = 10. 

The fact that we are computing the chromatic polynomial {v — —1) implies that the 
dimension of the transfer matrix should be equal to the number of non-crossing non-nearest- 
nei^'/i&or partitions (i.e., partitions in which no block contains a pair of nearest neighbors). 
We shall denote these numbers by TriRing(m), as we expect them to be equal to the dimen- 
sion of T for a cyclic triangular- lattice strip of width m (see Table 1). 

A further simplification comes from considering classes of partitions invariant under re- 
fiections with respect to the center of the strip. The number of such invariant classes is 
denoted by SqRing(m) in Table 1. Thus, we expect that these numbers are equal to the 
dimension of T for a cyclic square-lattice strip of width m. 

We have modified the perl code that we used for free/cylindrical boundary conditions 
[23, 24] to deal with the cychc boundary conditions. We have found that for m > 3, 
SqRing(m) is larger than the number of distinct eigenvalues SqRing'(m) found in [36, 37, 43] 
(see Table 1).^ In Section 2.3 below we show how the transfer matrix can be rewritten in a 
form that amounts to making each distinct eigenvalue appear with unit multiplicity; we also 
explain how the issue of degeneracies is linked to an underlying quantum-group symmetry. 

^Thc authors of Rcf. [43] use the notation NpXy.x to denote the number of distinct eigenvalues appearing 
in the chromatic polynomial for cyclic strips of square and triangular lattices of width Ly. 
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Chang and Shrock found a closed expression for the numbers SqRing' (m) using coloring- 
matrix methods [43] : 



LW2J , _ , X, 

sciRMg-w = 2(^-1)! E (2.15) 



k=0 



We notice that these numbers coincide with the sequence A025565 of [44], and we therefore 
have the alternative expression^ 

S,Ring'W^'"'E'(''";^)('",^;-'). (2.16) 



fe=0 



The asymptotic behavior for large m of the number of distinct eigenvalues was found by 
Chang and Shrock [43] to be 

SqRing'(m) - 3"'m-^/2 as m ^ oo (2.17) 



Remarks: 1) This transfer matrix is not invariant under the interchange of the top and 
bottom rows (i.e. i i'). The reason is that our transfer matrix only acts on the top row, 
leaving fixed the bottom one. Thus, we cannot reduce the size of the transfer matrix by 
including only connectivity states that are invariant under the transformation i ^ i'. 

2) With our transfer-matrix approach we have checked all previous results by Shrock and 
collaborators [36, 37, 39, 38]. 

3) Even though the final dimension of the transfer matrix T(m) is given by SqRing(m), 
the computations need to be performed in the larger space of non-crossing connectivities (of 
dimension C2m)- This fact limits the maximum width we are able to handle. 

m 

4) When we apply the operator Y[ \i' to a non-crossing connectivity v-p in (2.8), the 

i=l 

resulting connectivity state v-p' is not guaranteed to be non-crossing. However, as this step 
is the last one in the computation, we do not need to enlarge our vector space (to include 
all crossing connectivities). 

5) Wc have checked our results by comparing them to the appropriate chromatic poly- 
nomial with cylindrical boundary conditions 

where F and P denote free and periodic boundary conditions respectively. 

^The equivalence between (2.15) and (2.16) comes from the equality of their respective generating func- 
tions. In particular [43] , 

Vr?I-l = i:SqRing'(m)x™ 

m=l 
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2.2 Block structure of the transfer matrix 



The transfer matrix for a cyclic square-lattice strip of width m has a block structure. 
This property greatly simplifies the computation of its eigenvalues as we shall see below. 
Given a certain connectivity V of the top and bottom rows (sec, e.g., Figure 2), one can 
regard it as a top-row connectivity Vtop and a bottom-row connectivity ^bottom joined by £ 
bridges. Note that each top-row block is connected by at most one bridge to the bottom row 
and vice versa. Thus, i can take any integer value from to m. 

In order to obtain the matrix element T„ {, between two connectivity states w-p and w-p/, 
we should perform three operations 

1. Apply the vertical-bond operator V on the initial state w-p. 

2. Apply the horizontal-bond operator H, so we kill any nearest-neighbor connectivity 
state. 

3. Project the result onto the final state w-p>. 

Let us suppose we start from a connectivity state w-p with a certain bottom-row con- 
nectivity Pbottom, number of bridges i, and position of such bridges (all these characteristics 
modulo reflections with respect to the center of the strip). Then two important observations 
should be made: 



• The bottom-row connectivity of the initial state is not modified by the application of 
either matrix V (in step 1) or matrix H (in step 2). 

• The apphcation of the matrices V (in step 1) and H (in step 2) cannot increase the 
number of bridges i. Thus, we can only go to connectivity states V with the same 
or smaller number of bridges £' < £. Indeed, the relative position of the bridges is 
preserved. 

These two properties imply that the transfer matrix has a lower-triangular block form: 



/ T 



Tim) — 







Tm—l,m ^m— l,m— 1 



\ 



^0,0 / 



(2.19) 



(Note that we have arranged the connectivity states according to a decreasing number of 
bridges.) Furthermore, they also imply that each diagonal block T^^^ has a diagonal- block 
form 

/ Tf j A ... \ 



Ti^e. — 



'-i,e,i 

Te,e,2 ... 



V 











(2.20) 



where each sub-block T^,^j- is characterized by a certain bottom-row connectivity Pbottom and 
a position of the i bridges. Its dimension is given by the number of top-row connectivities 
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that are compatible with T'bottom, ^, and the relative positions of the £ bridges. The number 
of blocks A^^ (0 < ^ < m) for a square-lattice strip of width m is displayed in Table 2. 

In particular, this structure means that the characteristic polynomial of the full transfer 
matrix T can be factorized as follows 



X 



T) = n 



(2.21) 



1=0 

In practice, for each strip width 2 < m < 6, we have performed the following procedure: 

1. Using a perl script, we compute the full transfer matrix T(m) and the left u and 
right vectors Wid- These objects are obtained in the basis of non-crossing non-nearest- 
neighbor classes of connectivities that are invariant under reflections (with respect to 
the center of the strip). Thus, their dimension is given by SqRing(m). 

2. Given the full transfer matrix T(m), we compute the SqRing'(m) distinct eigenvalues 
Aj(g). To do this, wc split the transfer matrix into blocks T^/j-, each block charac- 
terized by a bottom-row connectivity and the number i and position of the bridges 
(modulo reflections). By diagonahzing these blocks we obtain the SqRing'(m) distinct 
eigenvalues Xi{q) and their multiplicities kiJ Indeed, 

SqRing'(m) 

ki ^ SqRing(m) (2.22) 

i=l 

3. The final goal is to express the chromatic polynomial in the form 

SqRing'(m) 

Pmxniq) - Yl <^ii^)^MT (2-23) 

i=l 

where the {Aj} arc the distinct eigenvalues of the transfer matrix T(m) and the {ctj} 
are some amplitudes wc have to determine. The calculation of the amplitudes can be 
achieved by solving SqRing'(m) linear equations of the type (2.23), where the l.h.s. 
(i.e., the true chromatic polynomials) have been obtained via Eq. (2.12) (i.e., by using 
the full transfer matrix and vectors). 

We have used this procedure to compute the transfer matrix up to width m = 6. The next 
case would involve the calculation of a symbolic matrix of dimension SqRing(m = 7) = 28940, 
which is unmanageable given our computer facilities. 

Remark. In the procedure described above one can obtain the eigenvalues \i{q) of the 
transfer matrix by considering the diagonal blocks T^^^ [c.f., (2.19)]. However, we still need 
the full transfer matrix in order to determine the amplitudes ai{q). 

^Unfortunately, we have not found any analytical formula for the multiplicities fcj. 
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2.3 Simplified computation of the transfer matrix 

The computations can be simplified, and taken to larger strip widths, by utilizing the 
following general property of the eigenvalue structure of T(m) : 

1. All distinct eigenvalues of the block Te^£ have the same amplitude, henceforth denoted 
a^^\ which is given by Eq. (2.25) below. 

This property is an analytical result originating from the representation theory of quantum 
groups [25, 26]. Before reviewing the arguments for it we shall pause to discuss its practical 
utility. 

Property 1 implies that the chromatic polynomial can be written as 

PmM = Y.^Y.{^ ' (2-24) 

^=0 1=1 

where is the number of distinct eigenvalues in the block 7)^^ (see Table 3). Note that we 
have labeled the eigenvalues as \^^\ with i being the corresponding number of bridges. 

One important consequence of Eq. (2.24) is that in order to obtain all the eigenvalues 
and amplitudes we do no longer need to compute the full transfer matrix T(m). Rather, 
it is sufficient to diagonalizc the blocks T^^i for £ = 0, 1, . . . , m to get the eigenvalues A^"'^ 
(j = 1, . . . , Lg) — assuming, of course, that we can determine the amplitudes a*^^^ by some 
other method. 

Property 1 and the exact expression of the amplitudes a*^^-* follow from the presence of 
a quantum-group symmetry in the square-lattice Potts model [25, 26]. As a first step, the 
Potts model is reformulated as a six- vertex model in a standard way [28], the cyclic boundary 
conditions being imposed by the insertion of an appropriate twist operator [26]. It is then 
realized that the vertex model transfer matrix Tyertex commutes with the generators of the 
quantum algebra UgSl(2) [25, 26], with deformation parameter q given by ^/q — q + 1/q. 
Decomposing the representation space into a direct sum over eigenspaces corresponding to 
each value of the vertex- model spin (which in turn is identified with i) gives Property 1, 
along with the expression 

= ( 4) ■ (2-25) 



^ 2 ^ 

where C/„(x) is the Chebyshev polynomial of the second kind, defined by [45], 



The form (2.25) has recently been rediscovered independently by Chang and Shrock [43], 
using arguments based on Temperley-Lieb algebra. Though the reasoning in [25, 26] and 
[43] may not constitute fully mathematically rigorous proofs, there is hardly any doubt that 
the results are correct. We have also verified Property 1 explicitly for strip widths m < 6. 
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The first coefficients a^^^ are given by 



= 1 (2.27a) 

a^^^ = q-l (2.27b) 

= q^-^q^i (2.27c) 

= g3 _ ^ _ ^ (2.27d) 

q;(4) = (g_ i)(g3 _g^2 ^g^_ (2.27e) 

Another useful characterization of the amphtudes is the following [43] : 



= n - (2-28) 



where Bn^ is the generalized Beraha number^ 

B^^ = 4cos2 = 2 + 2cos {^-^^ (2.29) 

and ^ ranges between and the strip width m. The standard Beraha numbers are S„ = Bn^ . 

In the case of the Potts model at general temperature, and generic values of g, the 
representation theory for Ug sl(2) also fixes the dimension of the spin £ eigenspaces as (^T^^) ~ 
(m-T-i) [25, 26]. This formula should be compared with Table 11 in [46]. Note however that 
this is not equal to the number of distinct eigenvalues Li in the spin £ sector of (2.24). 
There are two reasons for this: 1) we are here considering the special case of the chromatic 
polynomial, for which many eigenvalues vanish due to the non-nearest-neighbor constraint 
(see Table 3), and 2) by symmetrizing with respect to refiections of the strip we have paired 
eigenstates whose eigenvalues would otherwise have been degenerate. 

Let us stress that we have verified that Property 1 is in perfect agreement with all our 
exact transfer matrices up to strip width m — Actually, this empirical evidence for 
2 < m < 6 verifies a slightly stronger property, which we state here as a conjecture: 

2. Two distinct blocks T^^^ and T^/^^/ (with ^ ^ i'), cf. (2.19), have no common eigenvalue. 
In other words, there does not exist i, i', i, i' such that A^'^ = A^* ^ and £ ^ 

Note that only Property 1 is necessary for making our simplified computation of the transfer 
matrix work. However, we have rather strong reasons for believing that Property 2 is actually 
true in general, since it essentially states that the combined Ugsl(2) and reflection (in the 
case of the square lattice) symmetries provides a complete identification of all the symmetries 
in the problem. We leave the question whether Property 2 can be proved by quantum-group 
techniques aside for a separate publication. Note that Property 2 with Eq. (2.24) would 

m 

imply that — SqRing'(m), cf. Eqs. (2.15)-(2.16). And also note that Property 2 with 

Eq. (2.28) and the Beraha-Kahane- Weiss theorem [34] would imply that the only possible 

^Note that q = Bn ^ with n,j integer corresponds to the quantum-group deformation parameter q being 
a rational-order root of unity, q = e*'^-'/". 
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real or complex isolated limiting points for cyclic strips of the triangular and square lattices 
are odd-n generalized Beraha numbers.^ The fact that we observe exactly this structure of 
isolated limiting points also for m = 7 lends further credibility to Property 2. 

From the point of view of computing the transfer matrix, further efficiency is achieved 
by making a further empirical observation (based on the strips of widths 2 < m < 6): 

3. Let T^,e,jo be any sub-block of T^^, cf. (2.20), whose dimension is the largest possible 
(there may be more than one such sub- block) . Then the dimension of / jo exactly 
Li, and all of its eigenvalues are distinct. 

The dimensions of these sub-blocks are given by the in Table 2. 

The characterization of the sub-blocks T^/jg with maximal dimension is interesting from 

a computational point of view. In particular, we have found that, for a given value of the 
number of bridges i, the sub-blocks T^/jo split into two classes depending on whether the 
bottom-row connectivity ^bottom and the relative position of the £ bridges is symmetric or 
not under reflections. Let us denote those classes as Cgym and Cnon-sym, respectively. Then, 
we have empirically found for 2 < m < 8 that: a) All blocks belonging to class Cgym have 
the same dimension Lsym', b) All blocks belonging to class Cnon-sym have the same dimension 

-^'non-sym; and c) -/^non-sym 

> Lsym- Thus, any block belonging to the non-symmetric class 
Cnon-sym IS a maximal one with Lnon-sym = Li (see Table 3). The only exception to this rule 
occurs when the non-symmetric class is empty Cnon-sym ~ 0- In this particular case, any 
block belonging to Csym is a maximal one with = -^sym- For instance, for a strip of width 
m > 3, we can always choose as maximal blocks with k = 1, . . . ,m bridges those blocks 
characterized by Pbottom = 1 (i-e., no connections among the sites of the bottom rim) and 
bridges located at sites 1', . . . , k'. Finally, we can choose as maximal block with no bridges 
{£ — 0) the block with T^bottom = <^i',3'- 

Note that Property 3 is not very surprising at all when one recalls that jo labels the 
bottom-row connectivities. Since the transfer matrix acts on top-row connectivities only, 
one would not expect its spectrum to depend on jo (apart from the reflection parity effect 
observed above). To save space we refrain from turning this idea into a formal proof. 

Supposing Property 3 to be true in general, it suffices to diagonalize one maximal sub- 
block for each £ — 0, 1, . . . , m (e.g., those described above). This method has allowed us to 
extend our results up to m = 8. 

An important question is whether this procedure is able to give all the SqRing'(m) 
distinct eigenvalues of the transfer matrix T(m). For the strips of widths 2 < m < 6, we 
have exphcitly verified that this is the case, as we can handle the full transfer matrix (See 
Section 2.2). For widths m = 7, 8, the situation is apparently less clear. In Table 3 we have 
displayed the dimension Li{m) of the maximal block T^^ .,„(m) for a strip of width m with £ 
bridges. These numbers fully agree with the exact number of distinct eigenvalues np{m,£) 
in the corresponding ^-bridge subspace found by Chang and Shrock [43, Table 1]. Thus, we 
can guarantee that for m = 7, 8 we have not missed any transfer-matrix eigenvalue, and that 
our results are indeed exact. Furthermore, the total number of distinct eigenvalues we have 
found is equal to the predicted number SqFree'(m) [43]. 

^The position of the isolated limiting points depends on the boundary conditions of the strips. For square- 
lattice strips with free and cylindrical boundary conditions [22, 23] we have found that (at least up to widths 
m = 12) all Beraha numbers up to i?6 = 3 (and only those) can be isolated limiting points. Furthermore, 
we have found complex- conjugate pairs of complex isolated limiting points. 
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Remcirk. We stress that the results for widths 2 < m < 6 have been obtained completely 
ab initio, using the practical procedure described in Section 2.2. Thus, for these strips, we 
have not used any unproven property. For the strips of widths m = 7, 8, we have used the 
simplified method described above, which makes use of Properties 1 and 3. Property 2 is 
actually not needed and used in our practical procedure. 



3 Square-lattice chromatic polynomials with cyclic bound- 
ary conditions 

In this section we will analyze the results for the square-lattice strips of widths 2 < L < 8. 
Even though the two smaller strips (L = 2, 3) are well known in the literature, we will describe 
them in detail using our approach. The other cases published in the literature L — 4,5 will 
be briefly reviewed for completeness. Finally, our new results will be described in detail. A 
summary of the main characteristics of the limiting curves for these strips is displayed in 
Table 4. Finally, in Table 5 we list the real chromatic zeros for square- lattice strips of size 
Lf X {kL)p with 1 < A; < 10. 

For the larger L, we have not attempted to write the complete analytic expressions for the 
eigenvalues and amplitudes in this paper, as they are rather lengthy. The interested reader 
can find them in the MATHEMATICA file transf er4_sq.m available as part of the electronic 
version of this paper in the cond-mat archive. The analytic expressions for the sub-blocks 
Te,ijQ can be obtained by request from the authors. 



3.1 L = 2 

This case was first solved by Biggs, Damerell and Sands back in 1972 



Their solution 



IS 



^2xn(g) = {q' - 3g + 3)" + {q-l) [(3 - + (1 - qT] +q'-3q+l (3.1) 

Let us review our derivation of this result. The transfer matrix T(2) has dimension four. 
In the basis {5i,i'(^2,2', + ^2,2', ^1,2' + ^2,1', 1}, it takes the form 



T(2) 





1 





\ 




-1 


2-q 1 












1 2-q 







v 


1 


2{q-2) 2{q-2) 




-3? + 3 / 



(3.2) 



where we by vertical and horizontal lines show the block structure of this matrix (see below) . 
We have ordered the basis elements according to a decreasing number of bridges so that the 
matrix (3.2) shows clearly its lower-triangular block form. The right and left u vectors 
are given by 



= {A,2A,Q,A), w\ihA^q{q-l) 



w 



id 



(1,0,0,0) 



(3.3a) 
(3.3b) 



The transfer matrix (3.2) can be decomposed into three blocks: all of them are character- 
ized by the same bottom-row connectivity ('Pbottom = 1, i-e., there is no connection between 
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the sites 1' and 2'). In this case, the difference between blocks is solely due to the number 
of bridges £ connecting the top and bottom rows. 

The first block is one-dimensional and corresponds to the connectivity state 5i,i'52,2' 
(thus, £ = 2). The eigenvalue is 

Xe=2 = 1 (3.4) 

The third block is also one- dimensional and it corresponds to the singleton state 1. Thus, it 
is characterized by zero bridges {£ = 0). The eigenvalue is 

A,=o = g'-3g + 3 (3.5) 

Finally, the second block corresponds to the two remaining states, and it is characterized by 
a single bridge {£ = 1). The eigenvalues of this block are 

a£i = 3-g (3.6a) 
AfJi = 1-q (3.6b) 

The amplitudes associated to those eigenvalues are given by [c.f. (2.27)] 

ai=2 = a^'^^ (3.7a) 
Q;Ji\ = a(^) , J = 1, 2 (3.7b) 
ai=o = q;(°) = 1 (3.7c) 

where the a^'^ are given by (2.25), in full agreement with the exact result (3.1). As discussed 
in Section 2.3, the amplitude a^^^ is associated to every eigenvalue A^*^ coming from the block 

We have depicted the chromatic zeros for the square-lattice strips of width L = 2-p and 
lengths n = 10p,20p in Figure 3(a). We have also depicted the limiting curve B for this 
strip. ^° This curve was computed using the resultant method (see, e.g., [22]). It divides the 
complex ^-plane into four regions, and crosses the real g-axis at g = 0, 2. At this latter point 
it has a quadruple point. There are also two T points at g = 2 ± \/2. 

The isolated limiting points arc very easy to locate as the amplitudes have a simple 
analytic form (2.28): there is a single isolated limiting point at g = 1. 

3.2 L = ^ 

The solution for this strip was first found by Shrock and Tsai in 2000 [36]: 

^3xn(g) = (g^-5g2 + 6g-l)(-l)" + (g2-3g+l)[(g-l)"+(g-2)"+(g-4)T 



+(g-l) 



(-(g-2)T + 5]Ar 



j=6 



10 



E^" (3-8) 



i=9 



where the Aj for i = 9, 10 are the solutions of the equation 

- (g - 2)(g2 - 3g + b)x + (g - l)(g^ - 6g^ + 13g - 11) = (3.9) 

^"^The limiting curve for a square-lattice strip of width L = 2 and cyclic boundary conditions was first 
computed by Biggs, Damerell and Sands [35] . We include it here to compare it to the new results presented 
in this paper. 
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and those for i = 6, 7, 8 are the solutions of the third-order equation 

+ (2g2 - 9g + I2)x'^ + (g^ - lOg''^ + 36^^ - 56g + 2,l)x = 

(g-l)(g^-9g^ + 29g2-40g + 22) (3.10) 

We have re-derived the above result using our method. The transfer matrix is 20- 
dimensional, and can be split into eight blocks. The diagonalization of these blocks lead 
to the eigenvalues given in (3.8). We have also confirmed that the amplitudes are precisely 
those given by (3.8). 

There is a single block for £ = 3 bridges. Its eigenvalue is A^=3 — —1 and its amplitude 
is a^^\ We also find two blocks with i = 2. The corresponding eigenvalues are ~ 
g — 4, g — 2, g — 1, and their amplitude is a^'^\ There are three blocks with a single bridge 
{i = 1). The eigenvalues are X^l2i = ~{q ~ 2)^ and X^^S^''^''^^ are the solutions of the third- 
order equation (3.10). The amplitude is a^^^ for these four eigenvalues. Finally, there are 
two blocks with no bridges £ — 0; their eigenvalues are the solutions of the second-order 
equation (3.9), and a^^^ is their amplitude. 

The chromatic zeros for the square-lattice strips of width L = Sp and lengths n = 15p, 30p 
are displayed in Figure 3(b). The limiting curve B (computed with the resultant method) is 
also depicted in that figure. This curve was first obtained by Shrock and Tsai [36]. 

The hmiting curve crosses the real g-axis at g = 0, 2 and g fa 2.3365442725. We find 
ten T points at g « 1.7713445981 ± 1.7900810556 i, g ^ 2.0457165736 ± 1.6008334033 g ^ 
2.2701797354 ± 0.9704300675 i, q ^ 2.0401211769 ± 1.6256019349 i, and g ^ 2.3527863293 ± 
0.3787883673 i. We also find a single isolated limiting point at g = 1. 

3.3 L = 4 and L = 5 

The solution for L = 4 was obtained by Chang and Shrock [37]. We have re-derived 
their result by using our approach. The transfer matrix has dimension 94, and it can be 
split into 16 blocks. After obtaining the eigenvalues of all such blocks, we arrive at the 26 
distinct eigenvalues and amplitudes found by Chang and Shrock [37] . These eigenvalues and 
amplitudes have the structure described in Section 2.3. 

In Figure 3(c) we show the chromatic zeros for strips of sizes 4p x 20p and 4f x 40p, as well 
as the limiting curve B}^ The computation of this curve was done by means of the direct- 
search method (see e.g. [22]). Thus, the number of relevant features reported below should 
be regarded as lower bounds. The limiting curve B crosses the real g-axis at three points: g = 
0,2, andg 2.4928455591. We have found six T points at g 2.2192166698±1.5778851486i, 
g ^ 2.3408741083 ± 1.3290891386 i, and g ^ 2.5035997660 ± 0.4577795448 i. As in the 
previous strips, there is a single isolated limiting point at g = 1. 

The solution for L = 5 was found by Chang and Shrock [38]. The full transfer matrix 
has dimension 614, and it can be spht into 44 blocks. From them we have obtained the 
70 distinct eigenvalues and their associated amplitudes. Again we find the eigenvalue and 
amplitude structure discussed in Section 2.3. 

Although the limiting curve for this lattice was not included in the original reference by Chang and 
Shrock [37], it appeared in the subsequent paper [38] by the same authors. 
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In Figure 3(d) we depict the chromatic zeros for the strips of length n — 25p, SOp- We 
also plot the limiting curve B, which was computed using the direct-search method. This 
curve B crosses the real g-axis at g = 0,2 and q ^ 2.5823854661. We have found six 
T points at g fti 2.3986403427 ± 1.4510688289 ^, q ^ 2.5948107219 ± 0.4767225937i, and 
q 2.2911897201 ± 1.64934412051 Finally, there is a single isolated hmiting point at g = 1. 

3.4 L = 6 

The full transfer matrix T(6) has dimension 4028, and can be split into 105 blocks. After 
performing the diagonalization of these blocks, we have found 192 distinct eigenvalues (in 
agreement with Eq. (2.15)). The solution can be written in the form 

192 

Pexniq) = Y,a,{q)\i{qr (3.11) 

i=l 

We find several simple eigenvalues: e.g. A = 1, A = 4 — g, 3 — g, 2 — g, 1 — g. The 

other eigenvalues come from solving polynomic equations up to order 30. We find the same 
structure of the eigenvalues and amplitudes in terms of the number of bridges £ as for the 
other strips. 

In Figure 4(a) we show the chromatic zeros for the lengths 30p and 60p, and the limiting 
curve B. This curve has been computed using the direct-search method. Thus, the number 
of T points and the other features discussed below should be taken as lower bounds of the 
true values. 

The limiting curve B crosses the real g-axis at g = 0. 2 and g = go ~ 2.6460783059. 
We have found 10 T points at g ^ 2.3494422185 ± 1.6759458630 i, g ^ 2.3716213974 ± 
1.6415616784 i, g ^ 2.4271173045 ± 1.5496869588 i, g ^ 2.4988088784 ± 1.4046791200 i, and 
g « 2.6565980229 ± 0.4807158507 i. There are two small triangular-like regions delimited 
by the six T points around g ^ 2.37 ± 1.64i. In total, there are six enclosed regions in the 
complex g-plane; only two of them have non-void intersection with the real g-axis. 

There are two isolated limiting points at g = 1 and g = i?5, as the value of go is larger 
than B5. This is the smallest width for which we find two isolated limiting points among 
the square- lattice strips considered in this paper. 

We remark that the limiting curve crosses the real g-axis at go ~ 2.6460783059 > -B5. 
Thus, in principle there could be real zeros larger than g = (Sec Table 5).^'^ However, 
we computed the zeros of P6xn(?) up to n = 60 and did not find any real zeros larger 
than B5. This point is relevant in relation to a conjecture by Woodall [47], stating that no 
bipartite planar graph can have a real chromatic zero greater than g = 2. This conjecture 
was disproved in ref. [22], and modified in the following way: for any bipartite planar graph 
C, Poiq) > for aU real q > B^ [22, Conjecture 7.5]. 

^^This limiting curve was first obtained by Chang and Shrock [38]. 

^•^In refs. [22, 23] it was found that for square-lattice strips with cylindrical boundary conditions and even 
width L = 8, 10, 12, the limiting curve B does not cross the real g-axis. Indeed, for L = 4, 6, S crosses the 
real g-axis at qo < B^. No chromatic real root greater than B5 was found. 
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3.5 L = 7 



The full transfer matrix T(7) has dimension 28940, and can be spht into 294 blocks (see 
Table 2). The computation of such an enormous symbolic matrix is beyond our current 
computer facilities. However, using the simplified method described in Section 2.3, wc have 
been able to obtain the diagonal blocks T^/Jq (i-^-j ^^^^ each number of bridges £ = 0, . . . , 7, 
we choose among all possible sub-blocks of T^^^ one with the largest possible dimension given 
by in Table 3). 

The blocks T^,£jo themselves large: their dimensions range from 1 to 145. This fact 
severely limits the applicability of the method: we have been unable to compute the char- 
acteristic polynomial of any block of dimension larger than 100. Thus, with our current 
computer facilities, we cannot obtain closed formulas for the eigenvalues (as in the preceding 
cases). This implies that we cannot compute the chromatic polynomial for this strip family. 
However, we are able to compute the limiting curve (which depends only on the eigenval- 
ues provided none of the amplitudes au vanishes identically) by the direct-search method 
combined with numerical determination of the eigenvalues at any given q. 

In Figure 4(b) we show the zeros for a square-lattice strip of size 7f x 7p (obtained using 
Eq. (2.18) and previous results for cylindrical boundary conditions [22]). We also show the 
limiting curve B; it has been computed using the direct-search method. The computation 
of this curve is more involved than for the previous cases, as we do not have the chromatic 
zeros to guide our search. Thus, the number of T points and the other features discussed 
below should be taken as lower bounds of the true values. 

The limiting curve B crosses the real g-axis at g = 0, 2 and q = Qq ^ 2.6927280047. 
We have found 14 T points at q ^ 2.4008837606 ± 1.6805209516 i, q ^ 2.5993953677 ± 
1.2843636999 ?:, q ^ 2.7015691600 ± 0.4815451956 i, q ^ 2.7022195651 ± 0.9476253694 ?:, q ^ 
2.7126320142 ± 0.8591968848 i, q « 2.7152854276 ± 0.7936937430 i, and q f« 2.7248651978 ± 
0.845205995H. 

There are two tiny complex-conjugate features of the limiting curve that are not visible 
in Figure 4(b). Around q ^ 2.71 ± 0.85, we find two closed regions dehmited by the T points 
q ^ 2.7022195651 ± 0.9476253694 i and q ^ 2.7152854276 ± 0.7936937430 i. These regions 
are split in two by small lines connecting the T points q ~ 2.7248651978 ± 0.8452059951 i 
and q ~ 2.7126320142 ± 0.8591968848 i Thus, we find four small closed regions delimited 
by the latter four T points listed above. 

There are two isolated limiting points at g = 1 and q — B^, as the value of qo is larger 
than B^. 

3.6 L = 8 

The full transfer matrix T(8) has dimension 212370, and can be spht into 777 blocks 
(see Tabic 2). We have been unable to compute this huge symbolic matrix. By using the 
method described in Section 2.3, we have been able to obtain the nine diagonal blocks X^/jo- 
Several blocks T^/jo are still very large: their dimensions range from 1 to 385. Again, we 
cannot compute the symbolic characteristic polynomial of the largest sub-blocks, and hence, 
the chromatic polynomial. But the numerical computation of the eigenvalues is possible, but 
very time consuming. This is the reason we have not tried to determine the whole limiting 
curve B. The only exception is the determination of the value of go for this strip family: 
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qo 2.7287513178. Finally, we expect two isolated limiting points ai q — 1 and q — B^, as 
the value of qo is larger than B^. 

4 Transfer matrix for triangular- lattice strips with cyclic 
boundary conditions 

The results are very similar to those obtained for the square lattice in Section 2. The 
transfer matrix is written as 

T(m) = V(m) • H(m) (4.1) 

with the same formula for H(m) (2.4a), but a different one for the "vertical-bond" matrix 
V(m): 

V(m) = • Q^,^_i • Pm-i • . . . • Q2,i • Pi (4.2) 

In Figure 5 we show the construction of a 4 x 2 triangular-lattice strip. 

The chromatic polynomial is given by the same formulas as for the square lattice (2.5) / (2.12) . 
Again we choose to work with the connectivity basis w-p = H(m) ■ H'(m) ■ v-p, as w-p = 
for all partitions V that contain two nearest-neighbor sites (on cither the top or the bottom 
rows) in the same block. This basis has dimension TriRing(m) for a strip of width m (see 
Table 1). 

There is no reflection symmetry for a triangular-lattice strip, so the dimension of the 
transfer matrix is TriRing(m). However, the number of distinct eigenvalues should be the 
same as for the square lattice, viz. SqRing'(m) in Table 1, as the underlying Ug sl(2) quantum- 
group symmetry is the same [26]. That this is indeed the case has been checked explicitly 
by Chang and Shrock [43]. 

Also in the triangular-lattice case, the quantum-group symmetry ensures that the chro- 
matic polynomial for a strip of size m x n can be written as in (2.24). Again, this means 
that only odd-n generalized Beraha numbers can be isolated limiting points. 

The block structure of T(m) is similar to that of the square-lattice case (2.19) / (2.20) . The 
only difference is that the triangular lattice is not invariant under reflections with respect to 
the center of the strip. Thus, for a given number of bridges i, we expect that the number of 
blocks Nc is larger than for the square case (see Table 3). Indeed, for all m < 8, and for any 
£, we empirically find that A^^ coincides with the number of distinct eigenvalues Li coming 
from Tg^i, and that all the blocks have the same dimension (see Table 3). Thus, we 
conclude that for m < 8 

m 

= TriRing(m) (4.3) 

i=o 

^''For triangular-lattice strips with free and cylindrical boundary conditions, we have found that all Beraha 
numbers up to Bio are isolated limiting points [24] (at least up to widths m = 12). There are two different 
conjectures about how many Beraha numbers can be isolated limiting points for this lattice. Baxter's results 
[21]) imply that all Beraha numbers up to (and including) Bu should be isolated limiting points for large 
enough widths m. The reanalysis of Baxter's solution [24] led us to a new conjecture: only the Beraha 
numbers up to Bn (and possibly -612) can be isolated limiting points. In addition to these issues, wc have 
found that for strips with free and cylindrical there are no complex isolated limiting points up to widths 
m = 12. However, for "zig-zag" boundary conditions, we found complex isolated limiting points, and, more 
importantly, real isolated limiting points that are not Beraha numbers (e.g., q = 5/2). 
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and we conjecture that this equation holds for arbitrary values of m. 

Once again, the form (2.24) allows us to make computations for widths m > 6. Since all 
the sub-blocks T^/j- are expected to have the same dimension, we can choose any of them to 
perform the computations using the simplified method described in Section 2.3. 

Remarks: 1. As for the square- lattice case, we have checked our results with the help of 
the available results for triangular-lattice strips with cylindrical boundary conditions [c.f., 
Eq. (2.18)]. 

2. As for the square- lattice case, the transfer matrices for the strips of widths m < 6 have 
been obtained without using any unproven property. For m = 7, 8 we have used Properties 1 

and 3 of Section 2.3 to obtain all the eigenvalues. For each block Ti^£, we have obtained 
Li{m) distinct eigenvalues (See Table 3), where Li coincides with the predicted number of 
distinct eigenvalues for the sector with i bridges [41, Table 1]. 



5 Triangular-lattice chromatic polynomials with cyclic 
boundary conditions 

In this section we will analyze the results for the triangular-lattice strips of widths 2 < 
L < S. As in Section 3, we will describe our methods on the smaller widths L = 2, 3, 
briefly review (for completeness) the already- known cases (L = 4), and describe in detail the 
new results presented in this section {L > 5). In Table 6 we list the real chromatic zeros 
for triangular-lattice strips of size Lp x {kL)p with 1 < A; < 10. A summary of the main 
characteristics of the limiting curves for these strips is displayed in Table 7. 

The interested reader can flnd the analytic expressions for the eigenvalues and amplitudes 
for all lattices reported in this section in the mathematica flle transf er4_tri .m available as 
part of the electronic version of this paper in the cond-mat archive. The analytic expressions 
for the sub-blocks Ti^ejo can be obtained from the authors. 

5.1 L = 2 

This case was first solved by Tsai and Shrock [36] . The solution is 

^2xn(?) = (q' -Sq + l) + (q- + (g - 1) [A^ + A^] (5.1) 

where A3^4 are the solutions of the equation 

(g-2)2+ (2g-5)x + x2 = (5.2) 

Let us review our derivation of this result. The transfer matrix has dimension six. In the 
basis {5i,i'52,2', (^1,2', (^2,2', (^1,1', (^2,1', 1}, it takes the form 
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where the vertical and horizontal lines show the block structure of this matrix (see below). 
The right Wi^ and left u vectors are given by 



= {A,0,A,A,0,A), withA = g(g-l) (5.4a) 
wf^ = (1,0,0,0,0,0) (5.4b) 

Again, we have arranged the basis elements in decreasing order of number of bridges £, so 
the lower-triangular block form of T(2) becomes apparent. 

The full transfer matrix can be decomposed into four blocks: as in the square-lattice case 
with L — 2, all of them are characterized by the same bottom-row connectivity 'Pbottom = 
1, and differ only by the number of bridges i. The first block is one-dimensional and it 
corresponds to the connectivity state with i — 2 bridges. Its eigenvalue is 

Xe=2 = 1 (5.5) 

The fourth block is also one-dimensional and corresponds to £ = bridges. Its eigenvalue is 

A,=o = iq-2f (5.6) 

The second and third blocks are identical two-dimensional blocks corresponding to a single 
bridge £ = 1: for the former block it emerges at site 2', and for the latter one it emerges 
at site 1'. Their eigenvalues X^^li {j = 1)2) are given by the solutions of the second-order 
equation (5.2). Please note that the block Ti^i is block diagonal as discussed in Section 2.2. 

We have computed the amplitudes in the same way as for the square lattice. The results 
are given by [c.f. (2.27)] 

ai=2 = a^^^ (5.7a) 
aji^i = a^^\ i = l,2 (5.7b) 
a^^o = a^^^ = 1 (5.7c) 



These formulas agree with the already known exact results (5.1) / (5.2). Again, the amplitude 
associated to each eigenvalue A^"''' coming from block T^^i is a^^\ 

We have depicted the chromatic zeros for the triangular-lattice strips of width L — 2f 
and lengths n — 1 Op, 2 Op in Figure 6(a). The hmiting curve B for this strip is also shown. 
This curve has been obtained using the resultant method. It crosses the real g-axis at 
g = 0, 2, 3. At this latter point it has a quadruple point. The limiting curve divides the 
complex g-plane into three closed regions. Finally, there are two isolated limiting points at 
q — 1 and q — B^. 



5.2 L = S 

This case was first solved by Chang and Shrock [39], who also computed the limiting 
curve. We have re-derived their result using our method. The full transfer matrix has 
dimension 30. It can be split into ten different blocks. These blocks provide the ten distinct 

^^This limiting curve was first obtained by Shrock and Tsai [36]. 
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eigenvalues of the transfer matrix. Following a similar procedure as for the square-lattice 
strips we obtained the corresponding amplitudes. 

There is a single block with £ = 3 bridges: this is one-dimensional and it corresponds to 
the eigenvalue 

Xe=3 = -1 (5.8) 

and its amplitude is a'^^\ The block T2^2 contains three identical sub-blocks, each of them 
leading to three distinct eigenvalues X^~2''^'^^- The simplest one is 

= ? - 2 , (5.9) 

while the other two are the solutions of the second-order equation 

x'^ -x{2q-7)+q'^ -5q + 6 = 0. (5.10) 

The amplitude associated to these eigenvalues is a^'^\ We have also found four (identical) 
sub- blocks with a single bridge {i = 1). Each of them has four eigenvalues X^ll^- The simplest 
one is 

a£i = -q'^ + 5q-6 (5.11) 
and the other three are the solutions of the third-order equation 

= x^ + {2q^-12q + 19)x^ + {q^-llq^ + A6q^-85q + 5S)x 

-q^ + llg^ - 48g^ + lOV - 112g + 48 (5.12) 

The amplitude a^^^ corresponds to all of them. Finally, there are two (identical) sub-blocks 
characterized by no bridges {£ — 0). Their eigenvalues are the solutions of the second-order 
equation 

- a;(g^ - 7g2 + 18g - 17) + - + 30g2 - 44g + 24 = (5.13) 

and their amplitudes are a'-^-' = 1. 

In Figure 6(b) we show the chromatic zeros for the strips of sizes 3f x 15p and 3f x 15p, and 
the limiting curve B. This latter curve was obtained using the resultant method. It crosses 
the real g-axis at g = 0, 2, 3. At this latter point the curve has a quadruple point. There are 
eight T points at q ^ 2.2532978656 ± 1.6576481039 i, q ^ 3.1233767820 ± 1.1226454338 i, 
q « 3.1427903492 ± 0.3662742838 i, and q fti 3.0244195343 ± 0.1322645045 i. Finally we find 
two isolated limiting points at q — 1 and q — Br,. 

5.3 L = A 

This case was first solved by Chang and Shrock [39], who also computed the limiting 
curve. We have re-derived their result by using our general procedure. The full transfer 
matrix has dimension 178, and it can be split into 26 blocks. 

In Figure 6(c) we show the chromatic zeros for the strips of sizes 4f x 20p and 4f x 40p, 
and the limiting curve B}^ This latter curve was obtained using the direct-search method. 

^^The limiting curve for this strip was first computed by Chang and Shrock [39]. 
^''The hmiting curve for this strip was first computed by Chang and Shrock [39]. 
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Thus, the number of endpoints and other features reported below should be interpreted as 
lower bounds of the true values. 

The limiting curve B crosses the real g-axis at g = 0,2,3 and q ^ 3.2281262261. We 
have found eight T points at q ^ 3.3762587794 ± 0.9748801876 i, q « 3.3525265561 ± 
0.5771942145 i, q ^ 3.3542356970 ± 0.5945094569 i, and q ^ 1.4763063347 ± 2.03481318081 
There are two isolated limiting points at q — 1, B5. 

5.4 L = 5 

After we completed the computation of this case we learned that it had been indepen- 
dently solved by Chang and Shrock [41]. Our result based on our transfer-matrix method 
fully agrees with theirs. 

The full transfer matrix is 1158 dimensional, and can be split into 70 blocks. After 
the diagonalization of such blocks we find 70 different eigenvalues. Again, we find the 
eigenvalue /amplitude structure discussed in Section 2.3. 

In Figure 7(a) we show the chromatic zeros for the strips of sizes 5f x 25p and 5f x 50p, 
and the limiting curve B. This latter curve was computed via the direct-search method.^® 
The curve B has a more involved shape than for the other triangular-lattice strips with 
smaller width. It crosses the real g-axis at q = 0, 2, 3 and q ~ 3.3324469422. We have 
found 16 T points at q ^ 1.0114299869 ± 2.1132511328 i, q ft; 1.1643009748 ± 2.2002545351 
q ^ 2.5640974594 ± 2.0365055771 i, q ^ 2.5843109937 ± 1.9725597103 i, q w 2.6549519688 ± 
1.9693915215 i, q ^ 3.4850579943 ± 0.9796767236 i, q ^ 3.5231420699 ± 0.7322113419 i, and 
q ~ 3.5309588080 ± 0.8827974592 z. The unusual features arc a) two tiny bulb-like regions 
protruding from the T points at g pa 1.0114299869 ± 2.1132511328i; and b) two small 
triangular-shaped regions dehmited by the six T points around q ^ 2.60 ± 2.00 i 

Finally, there are three isolated hmiting points aiq — 1,-65 and B7. This is the triangular- 
lattice strip with smallest width for which we find an isolated limiting point at q — B^. 

5.5 L = 6 

The full transfer matrix is 7986 dimensional, and can be split into 192 blocks. The 
diagonalization of these blocks leads to the 192 different eigenvalues, in agreement with the 
formula (2.15). The simplest eigenvalue is A = 1 which comes from the block Tg g; thus, its 
amplitude is a^^\ The other eigenvalues are given by the solution of polynomic equations 
up to order 55. We find the block structure of eigenvalues and amphtudes described in 
Section 2.3. 

In Figure 7(b) we show the chromatic zeros for the strips of sizes 6f x 30p and 6f x 60p, 
and the limiting curve B. This latter curve was computed via the direct-search method. This 
curve crosses the real g-axis at g = 0, 2, 3 and q ~ 3.4062849655. We have found 10 T points 
at g fa 0.8028018425±2.2224309039i, q ^ 2.2195798331±2.3066154729 q ^ 3.6219823536± 
0.8757610205 i, q ^ 3.6246351106 ± 0.8660035657 i, and q ^ 3.6870064797 ± 0.8100610005 i. 
There are two bulb-like regions emerging from the last T points {q 3.6870064797 ± 

^^The computation of this limiting curve is new, as it was not published in the (first) preprint version of 
Ref. [41]. 
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0.8100610005 i). We have also found two endpoints at^^ q ^ 0.5716823545 ± 1.98961397731 
These points correspond to two complex-conjugate small branches protruding from the T 
points at g ^ 0.8028018425 ± 2.2224309039 i. This is first triangular-lattice strip for which 
we find endpoints. Finally, there are three isolated limiting points at g = 1, and Bt. 

5.6 L = 7 

The full transfer matrix T(7) has dimension 57346, and can be split into 534 blocks (see 
Table 3). By making use of the simplified method described in Section 2.3, we have obtained 
the relevant sub-blocks T^.^jo from which all distinct eigenvalues can be computed. 

Unfortunately, the sub- blocks T^/jo ^^^e rather large dimensions, ranging from 1 to 145. 
As for the square-lattice strip with L — 7, we have been unable to compute symbolically the 
characteristic polynomial of some of these blocks. Thus, the computation of the chromatic 
polynomial for finite strips of width L = 7 is beyond our computer facilities. However, as 
we can compute the eigenvalues and we know the amplitudes, we have been able to obtain 
the limiting curve. 

In Figure 7(c) wc show the zeros for a triangular-lattice strip of size 7f x 7p (obtained 
using Eq. (2.18) and previous results for cylindrical boundary conditions [24]). We also 
show the limiting curve B; it has been computed using the direct-search method. The same 
observations as for the square lattice with L — 7 hold in this case. 

The hmiting curve B crosses the real g-axis at q — 0, 2, 3, Bg and q — qo !v 3.4682618071. 
This is the first case for which wc find five different "phases" on the real g-axis. Wc have 
found 12 T points at q ^ 0.4951745217±2.1501328665i, q ^ 1.9229108754 ± 2.4842273215 i, 
q ^ 3.5014960312 ± 0.2613477419 q « 3.6121905140 ± 1.1057902185 i, q ^ 3.6937370146 ± 
0.8598201628 and q w 3.8081077493 ± 0.7447122915 i As for the L = 6 case, we find two 
complex-conjugate endpoints at g 0.3195792189 ± 1.8836441569.^° There are 3 isolated 
limiting points at q — 1, q — B^, and q — B^, as the value of go is smaller than Bg. 

5.7 L = 8 

The full transfer matrix T(8) has dimension 424206, and can be spht into 1500 blocks (see 
Table 3). We have computed the relevant nine blocks ^ from which all distinct eigenvalues 
can be obtained. As for L = 7, several of these blocks have large dimensions (up to 385). 
Thus, we have been unable to obtain the symbolic form of their characteristic polynomial. 
We can nevertheless compute the numerical values of the corresponding eigenvalues, and 
thus the limiting curve B. The computation of this curve is very time consuming, whence we 
have focused on obtaining the point go with great precision. Our result is go ~ 3.5134782609. 
As this value is greater than Bs, we can observe a "phase" characterized by £ = 4 bridges. 

""^^Unfortunately we have been unable to use the resultant method to verify the existence of these endpoints. 
Instead, we have used the direct-search method to locate them. The endpoint character of these points is 
suggested by the following observation. If Ai and A2 are two (equimodular) dominant eigenvalues, we can 
define the (real) parameter t = tan(6'/2), where 9 is given by A1/A2 = e'" (sec [22, Section 4.1.1]). Then, an 
endpoint is characterized by f = 0. At the points q w 0.5716823545 ± 1.9896139773 2, we have found that 
t « 1.8 X 10^*^. This computation strongly suggest that these points are indeed endpoints. 

^•^In this case, the parameter t = tan(0/2) (with 9 given by A1/A2 = e*^) takes the value t « 3.2 x 10~®. 
This small value suggests that these points are indeed endpoints. 
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Finally, we expect 3 isolated limiting points ai q — 1, q — Br,, and q — Bt, as the value 
of qo is smaller than Bg. 

6 Discussion 

We have empirically found an interesting structure of the eigenvalues and amplitudes of 
the transfer matrix in terms of the number of bridges of the connectivity states i. First, 
the transfer matrix has a lower-triangular block form (2.19). The blocks on the diagonal 
Ti^i in turn have a diagonal block form (2.20) in terms of the sub-blocks T^^ j. The sub- 
blocks Ti^ij are characterized by the bottom-row connectivity state and the number £ and 
position of the bridges (modulo reflection symmetry with respect to the center of the strip for 
the square-lattice case). Second, all the eigenvalues of the transfer matrix can be obtained 
by diagonalizing, for each i, only the largest of the sub-blocks. The final set of distinct 
eigenvalues can be then split into disjoint classes labeled solely by the number of bridges 
i (i.e., if a given eigenvalue appears in several sub-blocks, all of these sub-blocks belong to 
block Ti^e with i bridges). Finally, all the eigenvalues belonging to the class of i bridges have 
the common amplitude a^^K 

6.1 Phase diagram 

Motivated by this structure we have tried to give a topological meaning to the closed 
regions into which the limiting curve divides the complex g-plane. In other words, in each 
region of the g-plane we want to know to which block T^_^ the dominant eigenvalue belongs. 

We first discuss the regions having a non-empty intersection I with the real g-axis. We 
invariably find that: 

A. 7 is a simple interval in R* = R U {oo}, which is either 7„ = (-B2„, B2n+2) for n- — 
1,2,..., rimax - 1, where n„,^^ is determined by B2„^^^ < qo{L) < B2„^^^+2, or 1^^^ = 
iB2n^^, qo{L)), or Io = R*\ [0, qo{L)]. 

B. In the region containing 7„, with n — 0, 1, . . . ,ninax, the dominant eigenvalue comes 
from the block with i — n bridges. Its amplitude is thus a^^^. 

We have verified the Properties A-B above for all the cyclic strips of the square and 
triangular lattices considered in Sections 3 and 5. We conjecture that they hold true for 
any strip width L, for the appropriate value of qoiL), and in the thermodynamic limit with 
qo{L) replaced by go- [For a discussion on the values of go, see the Introduction. See also 
Figures 8-9 for plots of the limiting curves with 2 < L < 7.] 

In Ref. [26] a field theoretical mechanism was exhibited according to which the special 
role of the Beraha numbers in the case of free boundary conditions actually applies to the 
whole BK phase. We therefore finally conjecture that Properties A-B above hold within the 
entire BK phase, with qo{L) replaced by the boundary of the BK phase (in finite size, or in 
the 7/ — > oo limit). [For the extent of the BK phase in the thermodynamic limit, see the 
Introduction] . 

Property B should be compared with Saleur's study of the square-lattice Potts model 
along the curve v — —^/q, with < g < 4, i.e., the so-called unphysical self-dual fine. In this 
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case one has rimax = oo, and it can be proved that the dominant eigenvalue for q E In indeed 
occurs in the spin n sector [26] . By renormahzation group arguments one would expect this 
to extend to the entire BK phase. Another argument, asymptotically valid for L — cxo, is 
that for q E In the effective central charge is maximal in the n-bridge sector (see below). 

As to the regions inside the limiting curve that have an empty intersection with the 
real g-axis, we have found that in most cases (with 3 < L < 7) the dominant eigenvalue 
corresponds to the class i — 1. There are however some exceptions to this rule: 

• For the square lattice: For L — 7, the eigenvalues which are dominant inside the four 
tiny closed regions around q 2.71 ± 0.85 are characterized by either ^ = or £ = 2. 

• For the triangular lattice: For L — 5, the tiny bulb-like regions emerging from T points 
at g 1.0114299869 ± 2.1132511328 i, are characterized by £ = 2 (note that those 
regions are completely surrounded by an ^ = 1 phase). Furthermore, the two small 
triangular-shaped regions around q ^ 2.60 ± 2.00 i are characterized by = 0. For L = 
6, the two bulb-like regions emerging from T points at g ~ 3. 6870064797±0. 810061005 i 
are characterized by £ = 0. 

Therefore, we find that most of the limiting curves B corresponds to crossings (in modulus) 
of dominant eigenvalues coming from different blocks T^^^. However, there are some excep- 
tions to this general rule: we have found (in most cases, small) parts of B that correspond 
to crossings of dominant eigenvalues coming from the same block. For instance, the lim- 
iting curve for the square-lattice strip of width L = 2 contains a vertical segment which 
corresponds to crossings of two eigenvalues belonging to the i = 1 block. Smaller arcs with 
the same origin are also found for square- lattice strips with L < 3 (except L = 6), and for 
triangular- lattice strips with L > 5. 

6.2 Free energy and central charge 

In order to gain further insight into the phase diagram we have computed the free energy 
per spin fL{q) = 7; logA*(g) for real q for strips of width L and infinite length, where A* 
denotes the dominant eigenvalue of the transfer matrix. 

For L < 8 this was done using the symbolic technique described in Sections 2-5. Another 
technique enabled us to extend these results to the range L < 11. It consists in numerically 
diagonalizing the transfer matrix for a strip with fully free boundary conditions, but imposing 
the constraint that £ disjoint clusters propagate along the transfer direction. For each q, A*(g) 
was selected as the numerically largest among the L + 1 eigenvalues which are dominant with 
respect to each possible f-sector (^ = 0, 1, . . . , L). We have verified that for L < 8 the two 
techniques give identical results. 

In Figure 10(a) [resp. (b)] we plot /^(q') as a function of q for square-lattice (resp. 
triangular-lattice) strips with 2 < L < 11 (resp. 2 < L < 10). It is also useful to use 
instead of q, the (real) parameter n defined by 

q = 4cos (^^y , n>2 . (6.1) 

Special behavior (phase transitions, vanishing amplitudes, level crossings) is expected when 
n takes rational values. This is due to the special form (2.28) of the amplitudes a^^\ and to 
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the particularities of the representation theory of U^sl(2) when ^ is a root of unity [25, 26]. 
The vertical dashed lines in Figures 10-13 correspond to even Beraha numbers q — B^n- 

For a critical system, conformal field theory predicts the following finite-size scaling 
behavior with free transverse boundary conditions [48, 49]:^^ 

= Au,.(«) + ^ + ^ + o(L-) , (6.2) 

where /buik(Q') and /surf(?) are respectively the bulk and surface contributions to the free 
energy in the thermodynamic limit, c{q) is the central charge, and the geometrical factor 
G = \ (resp. G = 2/^/3) for the square (rcsp. triangular [50]) lattice. Neglecting all higher- 
order corrections on the right-hand side, one can obtain approximations to /buik(9), /surf(9), 
and c{q) from fits involving three system sizes: L — L^iin, -^min + 1, -^min + 2. However, we find 
that the o(L~^) terms are important and our estimates for c{q) have strong finite-size-scaling 
effects. 

In order to obtain accurate estimators for c(g) we need to include the first correction-to- 
scaling terms in (6.2). We have used the following improved Ansatz: 

f r \ f ( \ ^ /surf(g) , 7rG'c(g) A B 

The L^'^ term is the expected (non-universal) contribution of the non-singular part of the free 
energy due to the free transverse boundary conditions (see e.g., the discussion in Privman's 
review [51]). The last term is a non- universal correction predicted by conformal field 
theory. 

We have obtained approximations to /buik(?), /surf(?), and c{q) from fits involving five 
system sizes: L = L^in + k with < /c < 4. The best approximations for /buik(Q'), with 
Lmin = 7 (resp. Linin = 6) for the square (resp. triangular) lattice, are displayed as (black) 
solid circles in Figure 10. The corresponding estimates for c{q) arc shown in Figure 11, 
and for /surf(9) in Figure 12. Finally, for the triangular lattice we have also tried a four- 
parameter Ansatz, obtained from (6.3) by fixing /buik(?) to its analytically known value [21]. 
The resulting estimates for /surf(Q') and c{q) are shown in Figure 13. 

The estimates for the bulk free energy /buik(Q') shown in Figure 10 are rather smooth as 
a function of q, except close to go- In the triangular-lattice case (b), we have also shown the 
exact value of /buik(g) with cy/mtinca/ boundary conditions, as obtained by Baxter [21]. The 
very good agreement supports the idea that /buik(?) is independent of (reasonable) boundary 
conditions. However, for g ~ go this agreement is not longer good, due to large correction-to- 
scaling effects. In particular, for the triangular lattice our numerical results for 3.5 g 4 

^^Formula (6.2) is valid for two-dimensional systems only for manifolds with Euler number x = [49]. If 
X 7^ 0, then an additional term ~ logi appears [49]. 

^^Notc that (q) is the free energy per spin, whereas the standard CFT formulae require normalization per 
unit area. To redraw Figure 5 in an undistorted way (i.e., with each triangular plaquette being equilateral) 
the time slices must make an angle a = 27r/3 with the transfer direction. The projected width of an L-spin 
strip is thus L sina = L/G lattice spacings. 

^^There are indeed many more non-universal corrections terms due to irrelevant operators and some of 
them might be more relevant than the term included in the Ansatz (6.3). However, we expect that the 
term can mimic the effect of such higher-order terms. 
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do not shed light on the question of which of Baxter's eigenvalues is dominant where, and 
hence on the value of q'o(tri) (See [24] for a more comprehensive discussion on this issue). 

The results for c(g) make sense within the critical BK phase, i.e., for < g < goi where 
?o(sq) = 3 and go(tri) = 2 + \/3 (but see the Introduction for a discussion of the latter 
value). Estimates for both lattice types, and for different values of 2 < Lmin < 7, are 
displayed in Figure 11. The variation with respect to Lmin is a measure of the size of the 
residual corrections to scaling. 

It should be noted that the effective central charge Ccff (i.e., the one observed on the 
above plots) depends on both the "real" central charge Creai (i-e., the one computed from 
field theory) and the critical exponents corresponding to the insertion of I bridges. The 
latter are actually surface critical exponents, since the transverse boundary conditions are 
free. We have for the BK phase [52] 

C,^ = (6.4) 

n 

xe = (6.5) 

and since i — [n/2\ from Property A of Section 6.1, we find 

9zL~ 1 , (n-l)^-4(LtJ^-(n-l)LtJ) 
Ces{q) = c - 24:Xe = 1 - 6 ^ ^ . (6.6) 

This continuous oscillatory function is shown in Figures 11 and 13(b). Note that for n an 
even integer one has — 1 — 6/n, while for n an odd integer one has Cefr = 1 independently 
of n. However, exactly for integer n we do not expect (6.6) to hold (see below). Note that 
for free longitudinal boundary conditions by contrast, only the i = sector contributes, and 
thus the effective central charge (6.6) is equal to the BK value (6.4). 

Finally, for vahics of q outside the limiting curve we generally find that the effective 
central charge is close to zero, consistent with non-critical behavior. Indeed, with a finite 
bulk correlation length ^, the scahng form (6.3) must be replaced by 

Mq) = hMl) + + O (e-^/^ . (6.7) 

Note however that for the triangular lattice, wc expect the regime go(tri) < g < 4 to be 
critical [54, 52]. In particular the case q = 4 is equivalent to the three-coloring of the bonds 
of the hexagonal lattice [53], a critical model with c — 2 [55]. In addition, the 3-state zero- 
temperature square-lattice Potts antiferromagnetic model can be mapped to the F-model 
{alias the equal- weighted six- vertex model) [56, 57], which is a critical c = 1 theory. 

In the case of the triangular lattice, we have also tried a four-parameter Ansatz, obtained 
from (6.3) by fixing /buik(g) to its analytically known value [21]. In Figure 13(a) we have 
plotted the corresponding estimates for fsuriil), and in Figure 13(b) the estimates for the 
central charge. While we find a smooth behavior of /surf(?) for g < 3 and g ^ 4, there is 
a rather abrupt behavior in the regime 3 < g < 4. In that region we should expect strong 
finite-size-scaling corrections due to the proximity of several phase transitions. We conjecture 
that /surf(g) has a jump discontinuity at g = go, while /buik(g) is a continuous function of g 
(Baxter [21] explicitly showed that /buik(?) is continuous for the triangular lattice). 
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Our best numerical results for the central charge contained in Figures 11(a) (square 

lattice) and 13(b) (triangular lattice), show a remarkable agreement with the theoretical 
prediction for Cef[{q) (6.6) in the interval q G (0,2). We also find good agreement with the 
prediction c = in the non-critical region for q < —0.4, and q > 3.5 (resp. q > 4.5) for 
the square (resp. triangular) lattice. For q ^ and q qo there are strong corrections to 
scaling because the central charge is discontinuous at those points. Finally, we remark that 
at g = go we find excellent agreement with the theoretical predictions: c = 1.00(1) for the 
q = 3 model on the square lattice, and c = 2.01(5) for the q — 4 model on the triangular 
lattice.^"' 

6.3 Isolated limiting points and real crossings of the limiting curve 

We recall that isolated limiting points occur when the amplitude of the dominant eigen- 
value vanishes. In Tables 5 and 6 we show the real chromatic zeros for strips of size L-p x {kL)p 
with 1 < /c < 10. As explained in Ref. [22], the convergence to isolated limiting points is 
exponentially fast. 

First note that q;(°) = 1 never vanishes. Concerning the region containing /„ (with n — 
1,2, .. . ,ninax); wc remark that as q moves through the interval [B2n, B2n+2], ce'^"'^ increases 
monotonically from —1 to 1, passing through when q = -B2n+i- (This statement can be 
easily proved from the first of the definitions in (2.26).) In conjunction with the fact that 
all zeros of the Chebyshev polynomials of the second kind are real, the Properties A-B of 
Section 6.1 give the corollary: 

• An isolated limiting point must necessarily be real. 

• No isolated limiting point exists in the region containing Iq. For n = 1, 2, . . . , n^ax ~ 
1, the region containing has exactly one isolated limiting point which is -B2n+i- 
The region containing /nmax has the isolated limiting point -B2nmaxH-i provided that 

< QoiL). 

As before, we conjecture the validity of these properties for any L, and for the entire BK 
phase. 

By definition, the limiting curve crosses the real g-axis in qo{L). Further crossings are 
possible for < g < qo{L), provided that the two dominant eigenvalues become degenerate 
and that the corresponding amplitudes differ only by a sign. From (2.26) it is easy to see that 
the last condition is satisfied for q = with n = 1,2,..., nmax- That the first condition is 
also satisfied follows from the UgSl(2) symmetry [26], and we have also verified this directly 
for all strips studied in this paper. 

In the thermodynamic limit we have go(sq) = 3 and go(ti'i) = B12 = 2 + (as discussed 
in the Introduction). For the square lattice chromatic polynomial, we thus predict that the 
limiting curve crosses the real g-axis three times (at i?2, -B4 and Bq), dividing the real g-axis 
into three different phases (corresponding to £ = 0, 1, 2). Furthermore, there are two isolated 

"^^M q = qo we find that the dominant sector is £ = 0, thus there are no corrections due to the exponents 
xe (6.5). In addition, when performing the fits close to q = 4 for the triangular lattice (with the exact value 
of /buik('z)), we have found noticeable even-odd oscillations. The estimate c = 2.01(5) for g = 4 has been 
obtained by fitting data points with L = 4, 6, 8, 10. 



30 



limiting points (at and B5). For the triangular lattice, the limiting curve crosses the 
real g-axis seven times (at B2, B^. . . . , B12. and at -Boo)) dividing the real g-axis into seven 
different phases (corresponding to £ = 0, . . . , 6). Furthermore, there are five isolated limiting 
points (at -83,-65, • • -,-611). 



6.4 Cancellations for q = n integer 

One important point about the BK phase is that it is not defined for q = B^ with n 
integer [26]. In these cases there are cancellations among eigenvalues that give rise to a 
different physics. 

The most trivial case corresponds to q — B2 — Q. At this value the coefficients a^^'> (2.25) 
take the form 

„(% = 0)=^ ifmod(£.2) = j^^j 

1—1 otherwise 

We have seen that the eigenvalues become degenerate in pairs, one belonging to an odd-£ 
sector, and the other one to an even-£ sector. Thus there is a complete cancellation of all 
eigenvalues and the result is as expected, .^g„(0) = 0. 
At g = -B3 = 1, the coefficients take the form 

ff mod(£, 3) = 

aW(g = l)=<jo ffmod(£,3) = l (6.9) 

ff mod(£, 3) = 2 

The (non-zero) eigenvalues with non-zero amplitudes become degenerate forming 2n-tuples. 
Exactly n of them have a = 1, and the other half, a = —1. Thus, there is an exact 
cancellation of all contributions, giving Zc^iQ = 1) = 0. 
For q — B4 — 2, the coefficients take the form 

„(')(, = 2) = (l .fmod(Lf/2J.2) = j^j^j 

I — 1 otherwise 

In the triangular-lattice case we find a similar cancellation as for the q = 1 case: the non-zero 
eigenvalues with non-zero amplitude form 2n-tuples whose net contribution is zero. Thus, 
in this case we also have the expected result ^g„(? = 2) = 0. For the square lattice we find 
a non-critical behavior: there is an exact cancellation of all non-zero eigenvalues, except for 
those taking values A = ±1. The net contribution for a, L x n lattice is always of the form 
Zixn = 2 if n is even, and Zi^^n = if n is odd. 

The last trivial case corresponds to q = Bq = 3 for the triangular lattice. This is very 
similar to q = 2 for the square lattice: all non-zero eigenvalues cancel exactly, except those 
with A = ±1. Again, for each triangular lattice of size L x n the net contribution is given 
by Z^^n = 6 if mod(n, 3) = 0, and Z^^y^n = otherwise. 

For the square lattice the only non-trivial case is q = -85.^^ Up to -L = 5 we find that 
the dominant eigenvalue (belonging to the i = sector) actually contributes to the partition 

^^The value g = -Be = 3 is also non-trivial, as wc expect a critical point for the 3-statc antiferromagnet at 
zero temperature. However, we have found that the dominant eigenvalue (for all widths L < 8) belongs to 
the £ = sector, and a^°'> 7^ 0. 
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function. But for 6 < L < 8, we find that the dominant eigenvalue belongs to the 1 — 2 
sector, and a^'^\q = B^) — 0. The first non-zero contribution corresponds to a singlet 
belonging to A = 0. The numerical results are given in Table 8. 

For the triangular lattice, the value q = B^is also non-trivial. The dominant eigenvalue 
belongs to the 1 — 2 sector, but again cii^'^\q = B^) = 0. In this case, as L increases the first 
eigenvalue contributing to the partition function is deeper into the spectrum of the transfer 
matrix. For even L, this eigenvalue is complex and belongs io i. — 1, and for odd L, it is real 
and belongs to £ = 0. See Table 8 for the numerical results. 

For q = Bi, the dominant eigenvalue has i = for 2 < L < 4, and £ = 3 for 5 < L < 8. 
This latter case is the interesting one, as the corresponding amphtude vanishes. We find 
that the first eigenvalue that gives a net contribution to the partition function belongs to 
i = 0, and as L increases, its position in the transfer-matrix spectrum lowers rapidly. See 
Tabic 8 for the numerical results. 

These results show that exactly &t q = B^ the physics of the system is not described by 
the BK phase. In some cases g = 0, 1, 2 (and g = 3 for the triangular lattice) the partition 
function is trivial. In other cases, q — B^ and q — B^ (the latter one only for the triangular 
lattice), the free energy comes from a subdominant eigenvalue of the transfer matrix. We 
refrain from extracting the effective central charge from the results displayed in Table 8; as 
the data shows strong parity efi:ects and the number of data points is rather limited. 

A Alternative construction of the transfer matrix for 
square-lattice strips with cyclic boundary conditions 

We have already discussed in Section 2 how to build the transfer matrix for a square- 
lattice strip with cyclic boundary conditions: the idea was to leave the bottom row fixed 
and act only on the top row. This procedure leads to fairly large transfer matrices with a 
block structure. Each block is labeled by the bottom-row connectivity and the number and 
position of the bridges joining the top and bottom rows (modulo reflection symmetry) This 
transfer matrix is not invariant under the operation of interchanging the top and bottom 
rows. This latter operation can be viewed as a reflection with respect to the middle of the 
strip (in the longitudinal direction). For simplicity, we will call it a i-reflection. 

Smaller transfer matrices invariant under t-reflections can be built as follows. The idea 
is to obtain a transfer matrix acting on both the top and bottom rims. Thus, the matrices 
H and V take the form 

m— 1 m—1 

H(m) = W Qi,i+i ■ Qi',(i+i)' (A. la) 

i=l i=l 

m m 

V(m) = n^^-IIP^'- (^-Ib) 

i=l i=l 

These two matrices act on the top and bottom rows simultaneously, and are by construction 
t-invariant (i.e. invariant under a t-reflection i <-> i'). The t-invariant transfer matrix is 
deflned as 

T(m) = V(m) • H(m) . (A.2) 
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Note that as we are acting on both rims simultaneously, we can only obtain the partition 
function for strips of even length. Then, the partition function for a square-lattice strip of 
size m X 2n is given by 

P™x2n(g) = ■ H(m) ■ T(m)" • v,^ (A.3) 

where the matrices H and T arc given by (A.la)/(A.2). 

As we are considering the chromatic-polynomial case, the matrix H is a projector, and 
then it is more convenient to use the modified transfer matrix (see Section 2) 

V{m) = H(m) • V(m) • H(m) (A.4) 

and the connectivity-state basis 

wp = H{m) ■ v-p (A. 5) 

As in Section 2, we will consider hereafter this modified transfer matrix with the prime 
dropped. 

The fact the transfer matrix (A. 2) is now ^-invariant can be used to obtain a smaller 
transfer matrix. We can pass to a new basis consisting on connectivities that are either even 
or odd under i-reflection.^® In this new basis the transfer matrix takes the block-diagonal 
form 

- ( 'V f 1, ) (") 

where T_|_ (resp. T_) corresponds to the subspace which is even (resp. odd) under a t- 
reflection. The left and right vectors take the form 

< = Kd,+,0) (A.7a) 
= {ul,0) , (A.7b) 

as the start vector wid is clearly t-reflection invariant, and the action of on any t-reflection- 
odd connectivity is identically zero. Thus, the contribution of the t-reflection-odd subspace 
to the partition function vanishes, and wc can work entirely within the t-reflection-invariant 
subspace. The dimension of this subspace for a square-lattice of width m is denoted by 
SqRingT(m) in Table 9. We expect that for large enough width m, the dimension of the 
t-reflection-invariant subspace will be approximately one half of the original dimension. 

This transfer matrix also has a block structure. In this case, due to its invariance under 
t-reflection, the blocks are labeled merely by the number of bridges i joining the top and 
bottom rows. Thus, for a strip of width m, the transfer matrix T(m) can be split into 
i — m + 1 blocks, as in eq. (2.19): 



T(m) = 



/ fm,m ^0 ... \ 

T^m—l,m -^m— l,m— 1 • • • 



\ TQ,m-l ■ ■ ■ Tofi / 



(A.8) 



^^Let us denote {vi} the connectivity basis. Some of these elements are t-invariant; the rest can be 
grouped into pairs {va,Vb) that map into each other under a f-reflection. A basis for the even (i.e., t- 

rcflcction invariant) subspace is then given by the elements of the first set together with the combinations 
Va + Vb from the second set. A basis for the odd subspace is given by the combinations Va — Vb from the 
second set. 
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The difference with the method discussed in Section 2 is that now each diagonal block T^^i 
does not have in general any sub-block structure, like in 2.20. The dimension of the largest 
block of T is in general greater than the dimension of the largest block of T (see Tables 9 
and 2). As an example, for width m = 5, the largest block of the matrix T(5) has dimension 
21; but the dimension of the largest block of T(5) is 127. The fact that the dimension of the 
largest block for T grows very quickly with the strip width, makes the symbolic computation 
of its characteristic polynomial (and hence, of its eigenvalues) a difficult task. 

The simplest case is the strip of width m = 2. The basis (3.2) we used to construct the 
matrix T(2) is already i-reflection invariant. In this basis, the matrix T(2) takes the form 



7(2) 



with 



/ 1 





\ 


g-3 


g2 - 4g + 5 4 - 2g 







-1 


4 - 2g q^ -4:q + 5 













-3q + Sy J 


T42 


= 2(g-2)(g2-4g + 6) 







(A.9) 



(A.IO) 

The right Wid and left u vectors are equal to those given by (3.3). It is easy to verify that 



T(2) = T(2)^ 



(A.ll) 



This result is expected as both matrices are written in the same basis and each time we act 
with T(2) we add two layers to the strip. 

The first non-trivial case corresponds to the strip of width m = 3. In this case the t- 
reflection-invariant basis has dimension 15. We have found four blocks labeled by the number 
of bridges i: 



• £ — 3: This block is one-dimensional and the eigenvalue is A = 1. 

• i — 2: This block is four-dimensional and the eigenvalues are X — {q — 4)^, (g — l){q — 
4),(5-2)^and {q-iy. 

m i — 1: This block has dimension 10. One eigenvalue is X — {q — l){q^ — Qq^ + 13q — 11) 
and the other ones come from solving two different equations of third order. 

• i = 0: This block is three-dimensional. One eigenvalue is A = (g — 2)^ and the other 
two are the solutions of a second-order equation. 

Thus, we find 15 distinct eigenvalues. This number is larger than the number of distinct 
eigenvalues for this strip (see eq. (3.8)). Indeed, we have verified that the amplitudes cor- 
responding to five of the above eigenvalues are identically zero. As a matter of fact, the 
zero amplitudes correspond to the eigenvalues that are not positive definite for real q (e.g., 
A = {q — l){q — 4:) and {q — l){q^ — 6q'^ + 13q — ll)). This is expected as the transfer matrix T 
is physically equivalent to T^: hence, the former must contain the square of the eigenvalues 
of the latter. 

We have checked that this procedure works also for the square-lattice strip of width 
m — A and the result agrees with the one presented in the text. 
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Remark: This procedure can also be implemented for the triangular lattice. In this case, 
the relevant symmetry is a t-rcflcction followed by a (standard) reflection with respect to 
the center of the strip. This second step is needed to bring the diagonal bonds into the right 
position. 
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m 




TriRing 


SqRing 


SqRing' 
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4 
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3 
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30 


20 


10 


4 


1430 


178 


94 


26 


5 


16796 


1158 


614 


70 


6 


208012 


7986 


4028 


192 


7 


2674440 


57346 


28940 


534 


8 


35357670 


424206 


212370 


1500 


9 


477638700 


3210246 


1607246 


4246 


10 


6564120420 






12092 



Table 1: Dimensionality of the transfer matrix for cyclic boundary conditions. For each 

square-lattice width m we give the number of non-crossing partitions C2m, non-crossing 
non-nearest-neighbor partitions TriRing, the number of classes of partitions invariant under 
reflections SqRing, and the number of distinct eigenvalues SqRing' [43] . 
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8 


70 
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19 
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9 


179 


435 


528 


468 


318 


174 


72 


24 


5 1 


2204 



Table 2: Block structure of the transfer matrix of a cyclic square-lattice strip of width m as 
a function of number of bridges i. For each strip width m, we quote the number of blocks 
Nj( for a given number of bridges i, and the total number of blocks (# blocks). 
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835 
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334 


139 


43 


9 1 


4246 



Table 3: Block structure of the transfer matrix of cyclic square- and triangular-lattice strips 
of width m as a function of number of bridges £. For each strip width m, we quote the 
number of different eigenvalues L^ for a given number of bridges i, and the total number of 
distinct eigenvalues SqRing' (m). For the triangular lattice, the number of blocks Ni for a 
given value of i coincide with the number of distinct eigenvalues L^ for the same value of i. 





Eigenvalue-Crossing Curves B 
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Tabic 4: Summary of qualitative results for the square-lattice eigenvalue-crossing curves B 
and for the isolated limiting points of zeros. For each cyclic square-lattice strip considered 
in this paper, we give the number of connected components of B (# C), the number of 
endpoints (7^ E), the number of T points T), the number of double points D), and 
the number of enclosed regions (7^ ER). We also give the value go which is the largest real 
value where B intersects the real axis, and the number of real isolated limiting points of 
zeros RI). The symbol t indicates uncertain results. 
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4p X U)p 
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5p X 35p 


2 




5p X 40p 


2.000000000000 


2.566264491529 


5p X 45p 


2 




5p X 50p 


2.000000000000 


2.569485040710 


6p X 6p 


2.000004484676 


2.407498857052 


6p X 12p 


2.000000000000 


2.559303044172 


6f X 18p 


2.000000000000 


2.591819531465 


6p X 24p 


2.000000000000 


2.604283491675 


6p X 30p 


2.000000000000 


2.610364438319 


6p X 36p 


2.000000000000 


2.613696502926 


6p X 42p 


2.000000000000 


2.615612713769 


6p X 48p 


2.000000000000 


2.616718339556 


6p X 54p 


2.000000000000 


2.617340472059 


6p X 60p 


2.000000000000 


2.617677872418 


7p X 7p 


2 




Beraha 


2 


2.618033988750 



Table 5: Real zeros of the chromatic polynomial PmxniQ) of a square-lattice strip of width 
m and lengths n — k x m {k — 1, . . . , 10) with cyclic boundary conditions. A blank means 
that the zero in question is absent. The first two real zeros g = 0, 1 are exact on all lattices. 
For m = 2 we have only found these two first zeros. "Beraha" indicates the Beraha numbers 
£4 = 2 and B5 = (3 + y5)/2. 
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Lattice 


4tn Ziero 


otn Ziero 


otn Ziero 


/tn Zero 


2f 


X 2p 


3 








2f 


X 4p 


2.515844688131 


3 


3.408771911807 




2f 


X 6p 


2.640265063616 








2f 


X 8p 


2.616262917475 


2.751508037623 


3 




2p 


X lOp 


2.614657987429 


3 


3.138323546170 




2f 


X 12p 


2.621439264722 








2f 


X Hp 


2.616560224418 


2.887655738411 


3 




2f 


X 16p 


2.618700179921 


3 


3.083264711515 




2f 


X 18p 


2.617797220618 








2f 


X 20p 


2.618107109614 


2.926805351303 


3 




3f 


X 3p 


2.546602348484 








3f 


X 6p 


2.617993627116 








3f 


X 9p 


2.618034110186 


2.937234729089 






3f 


X 12p 


2.618033988011 








3f 


X 15p 


2.618033988754 


2.965427896259 






3f 


X 18p 


2.618033988750 








3f 


X 21p 


2.618033988750 


2.976132937129 






3f 


X 24p 


2.618033988750 








3f 


X 27p 


2.618033988750 


2.981771802850 






3f 


X ;-!0p 


2.()1803398875() 








4f 


X 4p 


2.617986010522 


3 


3.465246100723 




4f 


X Bp 


2.618033988761 


3 


3.230317951180 




4f 


X 12p 


2.618033988750 


3.000090191719 


3.165028046404 




4f 


X 16p 


2.618033988750 


3 


3.231342732991 




4f 


X 20p 


2.618033988750 


3 


3.219343048602 




4f 


X 24p 


2.618033988750 


3.000000003255 


3.200001385558 




4f 


X 28p 


2.618033988750 


3 


3.219217913659 




4f 


X 32p 


2.618033988750 


3 


3.217140972107 




4f 


X 36p 


2.618033988750 


3.000000000000 


3.210324007476 




4f 


X 40p 


2.618033988750 


3 


3.218488888204 




5f 


X 5p 


2.618033990394 


3 






5f 


X lOp 


2.618033988750 


3 


3.246585484861 




5f 


X 15p 


2.618033988750 


2.999999999573 


3.250186825428 


3.287751313682 


5f 


X 20p 


2.618033988750 


3 


3.246967522201 




5f 


X 25p 


2.618033988750 


3 


3.246982088945 


3.332102289607 


5f 


X 30p 


2.618033988750 


3.000000000000 


3.246977286183 




5f 


X 35p 


2.618033988750 


3 


3.246979668997 


3.329353767785 


5f 


X 40p 


2.618033988750 


3 


3.246979593498 




5f 


X 45p 


2.618033988750 


3.000000000000 


3.246979607399 


3.327447532663 


5f 


X 50p 


2.618033988750 


3 


3.246979603508 




6f 


X 6p 


2.618033988750 


3.001033705947 


3.125892136302 




6f 


X 12p 


2.618033988750 


3.000000000036 


3.246874398154 




6f 


X 18p 


z.OloUooyooYOU 


o . UUuUuUUUUUuU 






6f 


X 24p 




3 nonnnnnnnonn 






6f 


X 30p 


2.618033988750 


3.000000000000 


3.246979603717 




6f 


X 36p 


2.618033988750 


3.000000000000 


3.246979603717 




6f 


X 42p 


2.618033988750 


3.000000000000 


3.246979603717 




6f 


X 48p 


2.618033988750 


3.000000000000 


3.246979603717 




6f 


X 54p 


2.618033988750 


3.000000000000 


3.246979603717 




6f 


X 60p 


2.618033988750 


3.000000000000 


3.246979603717 




7f X 7p 


2.618033988750 


3 


3.247001348628 


3.404690481534 


Beraha 


2.618033988750 


3 


3.246979603717 


3.414213562373 



Table 6: Real zeros of the chromatic polynomial -Pmxn(?) of a triangular-lattice strip of width 
m and lengths n — k x m {k = 1, . . . , 10) with cyclic boundary conditions. A blank means 
that the zero in question is absent. The first three real zeros q — 0, 1, 2 are exact on all 
lattices. "Beraha" indicates the Beraha numbers B5 — {3 + -\/5)/2, Bq — 3, S7, and Bs- 
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Eigenvalue-Crossing Curves B 


Isolated romts 


Lattice 


#C 


#E 


#T 


#D 


#ER 


Qo 


#RI 


2r 


1 








1 


2 


3 


2 


3r 


1 





8 


1 


6 


3 


2 


4r 


It 


ot 


8t 


Ot 


5t 


3.2281262261 


2 


5r 


It 


ot 


16t 


Ot 


9t 


3.3324469422 


3 


6r 


It 


2t 


lot 


Ot 


5t 


3.4062849655 


3 


7f 


It 


2t 


12t 


ot 


6t 


3.4682618071 


3 


8f 












3.5134782609 


3 



Table 7: Summary of qualitative results for the triangular-lattice eigenvalue-crossing curves B 
and for the isolated limiting points of zeros. For each cyclic triangular-lattice strip considered 
in this paper, we give the number of connected components of B C), the number of 
endpoints (# E), the number of T points (# T), the number of double points (# D), and 
the number of enclosed regions ER). We also give the value go which is the largest real 
value where B intersects the real axis, and the number of real isolated limiting points of 
zeros RI). The symbol t indicates uncertain results. 



Q 


Lattice 


L 


^ / 




Square 


6 





0.28823784 






7 

8 







0.28514284 

0.28291082 




Triangular 


2 


1 


-0.24060591 






3 





-0.35114882 






4 


1 


-0.46727989 






5 





-0.48121183 






6 


1 


-0.48248246 






7 

8 




1 


-0.52290292 
-0.53899001 




Triangular 


5 





0.11980838 






6 





0.11458681 






7 





0.11124861 






8 





0.10886700 



Table 8: Real free energy for q — Bn with integer n. For several values oi q — Bn we list 

the free energy associated to the most relevant eigenvalue of the transfer matrix that gives 
a non-zero contribution to the partition function. For each value of q, the type and width L 
of the lattice strip, we list the sector i the eigenvalue belongs to and the free energy /. 
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m 


SqRing 


SqRingT 


# blocks 


max dim 


SqRing' 


2 


4 


4 


3 


2 


4 


3 


20 


15 


4 


7 


10 


4 


94 


60 


5 


25 


26 


5 


614 


342 


6 


127 


70 



Table 9: Block structure of the square-lattice transfer matrix acting on the top and bottom 
rows T (A.2). For each strip width m, we quote the dimension of the standard connectivity- 
state basis {w} SqRing and the dimension of the f-reflection-invariant basis SqRingT. We 
also give the number of blocks blocks) the transfer matrix T can be decomposed, the 
dimension of the largest block (max dim), and the number of distinct eigenvalues we find 
SqRing'. 
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Figure 1: Square lattice with cyclic boundary conditions of size 4x2. This lattice is obtained 
from a square lattice with free boundary conditions of size 4 x 3 by identifying the top and 
bottom rows i <-> i' {i — 1, ... ,4 — m). 



12 3 4 



^ i ^ • 

1' 2' 3' 4' 

Figure 2: Connectivity state for a cyclic strip of width 4. As an example, we show the 
connectivity state V — 5i',3'(52,4,4'- This state can be seen as a bottom-row connectivity state 
^bottom = Si',3' (solid line), a top-row connectivity state Ptop — <^2,4 (solid line), and a bridge 
connecting both rows 54,4/ (dashed line). 
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Zeros sq lattice = 2^ 



Zeros sq lattice = 3^ 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




-10 12 3 -10 12 3 

Re(q) Re(q) 

(c) (d) 



Figure 3: Limiting curves for square-lattice strips of width (a) L = 2, (b) L = 3, (c) L = 4, 
and (d) L = 5 with cychc boundary conditions. We also show the zeros for the strips 
Lp X (5L)p (black □) and Lp x (lOL)p (red o) for the same values of L. 
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Zeros sq lattice = 6jj Zeros sq lattice = 7jj 




-10 12 3 -10 12 3 

Re(q) Re(q) 

(a) (b) 



Figure 4: Limiting curves for square-lattice strips of width (a) L = 6, and (b) L = 7 with 
cychc boundary conditions. In (a) we show the zeros for the strips 6f x 30p (black □) and 
6f X 60p (red o). In (b) we show the chromatic zeros for the strip 7f x 7p (black □). 




Figure 5: Triangular lattice with cyclic boundary conditions of size 4x2. This lattice is 
obtained from a triangular lattice with free boundary conditions of size 4 x 3 by identifying 
the top and bottom rows z z' (z = 1, . . . , 4). 
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Zeros tri lattice = 2^ Zeros tri lattice = 3^ 




Re(q) Re(q) 

(a) (b) 

Zeros tri lattice L = 4„ 

X t\ 

I 1 ' 1 ' 1 — I 




I I I I I I I 

2 4 

Re(q) 

(c) 

Figure 6: Limiting curves for triangular-lattice strips of width (a) L = 2, (b) L = 3, and (c) 
L = 4 with cyclic boundary conditions. We also show the zeros for the strips Lp x (5L)p 
(black □) and Lp x (lOL)p (red o) for the same values of L. 
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Zeros tri lattice L - 5„ Zeros tri lattice L - 6„ 

X JX X JX 




Re(q) Re(q) 

(a) (b) 

Zeros tri lattice L = 7„ 

X t\ 

I 1 ' 1 ' 1 — I 




I I I I I I I 

2 4 

Re(q) 

(c) 

Figure 7: Limiting curves for triangular-lattice strips of width (a) L = 5, (b) L = 6, and (c) 
L = 7 with cyclic boundary conditions. We also show the zeros for the strips Lp x (5L)p 
(black □) and Lp x (lOL)p (red o) for L = 4,5. For L = 7 we show the chromatic zeros for 
the strip 7f x 7p (black □). 
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Limiting Curves Sq Lattice 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I I I I I I I I I I I I I L 



-10 12 3 

Re(q) 

Figure 8: Limiting curves for square- lattice strips of width L and cyclic boundary conditions. 
We show the results for L = 2 (black), L = 3 (red), L = 4 (green), L = 5 (blue), L = 6 
(light blue), and L = 7 (pink). The outermost (black) curve shows the limiting curve for a 
square- lattice strip of width L = 10 with cylindrical boundary conditions [23]. 
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Limiting Curves Tri Lattice 

\ ' \ ' r 




I I I I I I I 

2 4 

Re(q) 

Figure 9: Limiting curves for triangular-lattice strips of width L and cyclic boundary condi- 
tions. We show the results for L = 2 (black), L = 3 (red), L = 4 (green), L = 5 (blue), and 
L = 6 (light blue), and L = 7 (pink). The outermost (black) curve shows the infinite- volume 
limiting curve obtained by Baxter [21]. See [24] for a critical discussion on Baxter's limiting 
curve. 
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AF sq-lattice Potts 



AF tri-lattice Potts 




Figure 10: Estimates of the free energy for the square- (a) and triangular-lattice (b) Potts 
model for real q in the range — 1 < (? < 5. In each plot we show the values of the free energy 
obtained from the leading eigenvalue of the transfer matrix for strip widths L = 2 (black 
A), 3 (red o), 4 (blue x), 5 (pink +), 6 (orange •), 7 (green ■), 8 (violet A), 9 (red ♦), 
10 (blue *), and 11 (pink ©). We also show the extrapolated value of the free energy in 
the thermodynamic limit (black •) obtained by performing the fit (6.3) for Lmin = 7 (resp. 
Lmin = 6) for the square (resp. triangular) lattice. In (b) the solid (red) line shows the exact 
value of the triangular-lattice free energy found by Baxter [21]. The vertical (pink) dashed 
lines show even Beraha numbers q = = B^^^ (2.29). In (a) we show 1 < n < 3, and in 
(b), 1 < n < 4 and the limit oo (i.e., g = 4). 
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Figure 11: Estimates of the central charge for the zero-temperature Potts antiferromagnet 
on the square (a) and triangular (b) lattices for real q in the range — 1 < < 5. In each 
plot we show the values of the central charge obtained from the fits to the Ansatz (6.3) for 
-^min = 2 (black A), 3 (red o), 4 (blue x), 5 (pink +), 6 (orange •), and 7 (green ■). The 
solid (red) line shows the theoretical prediction for the effective central charge Ceg (6.6). The 
vertical (pink) dashed lines are as in Figure 10. 
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AF sq-lattice Potts AF tri-lattice Potts 
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q q 



(b) 



Figure 12: Estimates of the surface free energy /surf(Q') for the zero-temperature Potts anti- 
ferro magnet on the square (a) and triangular (b) lattices for real q in the range — 1 < q' < 5. 
In each plot we show the values of the surface free energy obtained from the fits to the 
Ansatz (6.3) for Lmm = 2 (black A), 3 (red o), 4 (blue x), 5 (pink +), 6 (orange •), and 7 
(green ■). The vertical (pink) dashed lines are as in Figure 10. 
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Figure 13: Estimates of the surface free energy /surf(Q') (a) and the central charge c(g) (b) 
on the triangular lattice when /buik(Q') in the Ansatz (6.3) is fixed to its analytically known 
value [21]. We show the estimates of /surf(9) and c(g) for real q in the range — 1 < g < 5 
and different values of L^in = 2 (black A), 3 (red o), 4 (blue x), 5 (pink +), 6 (orange •), 
and 7 (green ■). The solid (red) line in (b) shows the theoretical prediction for the effective 
central charge Ces (6.6). The vertical (pink) dashed lines are as in Figure 10. 
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